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Abstract 



The COSY-11 collaboration measured the pp pprj and pp — > ppr]' reactions in 
order to perform comparative studies of the interactions within the proton-proton-meson 
system. 

This thesis presents in detail the analysis of the pp — > pp-q' reaction which was measured 
at the proton beam momentum of 3.260 GeV/c. 

The elaboration results in differential distributions of squared invariant proton-proton {spp) 
and proton- r/' (spr;') masses, as well as in angular distributions and the total cross section 
at an excess energy of 16.4 MeV. 

The differential distributions Spp and Spr,' are compared to theoretical predictions and to 
the analogous spectra determined for the pp pprj reaction. 

The comparison of the results for the rj and rj' meson production rather excludes the hy- 
pothesis that the enhancement observed in the invariant mass distributions is due to the 
meson-proton interaction. 

Further, the shapes of the distributions do not favour any of the postulated theoretical 
models. 



Streszczenie 

W ramach grupy badawczej COSY-11 wykonano pomiary reakcji pp — > pprj oraz 
pp ppr]' w celu przeprowadzenia studiow porownawczych oddzialywania w ukladzie 
proton-proton-mezon. 

W rozprawie doktorskiej zaprezentowano analiz^ danych z pomiaru reakcji pp pprj' 
wykonanego z wykorzystaniem wi^zki protonowej o p^dzie 3.260 GeV/c. 
Wynikiem analizy przedstawionej w niniejszej pracy s^ rozniczkowe przekroje czynne w 
funkcji mas niezmienniczych proton-proton {spp) i proton- r/' (sp^'), rozklady k^towe oraz 
calkowity przekroj czynny dla energii wzbudzenia Q = 16.4 MeV. 

Rozklady rozniczkowe Spp oraz Sprf zostaly porownane z przewidywaniami teoretycznymi 
oraz analogicznymi widmami otrzymanymi dla reakcji pp pprj. 

Porownanie rezultatow otrzymanych dla produkcji mezonow rj i rj' pozwala na wykluczenie 
hipotezy, ze wzmocnienie obserwowane w widmach mas niezmienniczych jest powodowane 
oddzialywaniem mezonu z protonem. 

Co wi^cej, na podstawie ksztaltu otrzymanych dystrybucji nie da si^ rozstrzygn^c pomi^dzy 
poprawnosci^ zalozeh w postulowanych modelach teoretycznych. 
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1. Introduction 



The knowledge of the meson-nucleon interaction at the hadronic level is one of the ma- 
jor goals in nuclear physics, nowadays. Also, studies of meson structure and production 
mechanisms constitute a huge interest of nuclear and particle physicists. During the last 
decades many measurements produced interesting results [1], but various questions are still 
open. 

In the SU(3)-flavour scheme the rj and ?/ mesons belong to the nonet meson family of 
pseudoscalar mesons. The rj and ?/ mesons constitute a mixture of states: singlet-r/i and 
octet-r/g. 

It is important to stress that the strength of the prj and prj' interaction depends on the 
structure of the r] and rj' mesons and is directly related to the singlet-r^i and octet-r/g con- 
tributions in the wave functions of these mesons [2-4]. 

Taking into account the mixing angle the resulting contribution of the various quark 
flavours in the t] and 7]' wave function is almost the same. However, in spite of the postu- 
lated similar quark structure the r] and rj' mesons owe unexpectedly different features. 
The most drastic ones are: 

• rj and t]' mesons possess different masses, the 77(547) meson is almost two times lighter 
than the r/'(958) meson [5]. 

• The branching ratios of B and D<j mesons for the decays into the r]' meson are higher 
than for the decays into the rj meson, and deviate strongly from model predictions 
[6,7]. 

• There is no experimental evidence for baryonic resonance which decays by the emis- 
sion of an r]' meson [5,8] while e.g. the resonances N(1535) and N(1650) decay via 
the emission of the r/ meson with a significant probability [5]. 

These different features could imply a possible difference in the interaction of rj and r]' 
mesons with elementary particles, and indicate that also the production mechanism of 
both mesons in elementary particle collisions might vary. 

Due to the short life-time of the flavour neutral pseudoscalar mesons, experiments with 
meson beams or targets are difficult or hardly feasible. Therefore, the nucleon-meson in- 
teraction can be studied only via its influence on the cross sections of the reactions during 
which they are produced (e.g. NN NN Meson) [1]. 

Quantitative information about the interaction can be gained from the shape of the exci- 
tation functions for the pp — > pprj and pp — > pprj' reactions as well as from a comparison 
of those to the ppn^ system. 
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Up to now, only the proton interaction with pions and the r/-meson [1] was studied 
more exhaustively due to much more higher total cross sections in comparison to the cross 
section for the pp — > pprj' reaction. 

Besides the excitation function also differential distributions of invariant proton-proton 
and proton-meson masses constitute a sensitive tool for studies of the interaction within 
the meson-nucleon system. The distribution of the proton-77 invariant mass showed a 
clear enhancement in the region of small proton-?] relative momenta [9]. The observed 
effect could be explained by the non negligible role of the proton-ry interaction in the final 
state [10, 11], the admixture of higher partial waves during the r/ production [12], or the 
energy dependence of the production amplitude [13]. Using only the pp — > ppr] data, it is 
however not possible to justify or falsify any of the above mentioned hypothesis. 

The endeavour to explain the observed enhancement motivated the experiment and the 
analysis of the pp — > ppr]' reaction which is presented in this thesis. This experiment was 
performed using the COSY-11 detector setup [14-16], installed at the cooler synchrotron 
COSY [17] at the Research Centre Jiilich in Germany. It was conducted in an energy range 
close to the kinematical threshold for the t]' meson production, where the relative velocities 
of the produced particles are small. 

The analysis and results of the pp — > ppr]' reaction measurement, conducted in order to 
determine the distribution of events over the phase space for an excess energy range equal 
to the one measured before for the pp — > ppr] reaction are presented in this thesis. The 
measurement was performed at the nominal beam momentum of 3.257 GeV/c correspond- 
ing to the nominal excess energy for the ppr]' reaction equal to 15.5 MeV. The analysis of 
the ppr]' data was performed in a similar way as it has been done for the ppi] system. 
The comparison of the differential distributions for the proton-proton and for the proton- 
meson invariant masses in the rj and r]' production could help to judge about the validity 
of postulated theories concerning the observed enhancement and allows for a quantitative 
estimation of the relative strength of the proton-77 the proton-ry' interactions, provided 
that the effect is caused by the proton-meson interaction. 

In the following chapter the current status of the pp — > ppr] and pp — > ppr]' total cross 
section measurements is presented. Furthermore, the possibility of proton-r/ and proton-r]' 
interaction studies is discussed and the description of the observables used in the analysis 
described in this thesis is given. The results achieved for the pp ppr] measurement as 
well as the available theoretical descriptions of the results are presented. 
Chapter [3] is devoted to the presentation of the experimental facility used to perform the 
pp — > ppr]' measurement. In the next chapter, the methods used for the calibration of the 
COSY-11 detectors and their relative geometrical settings are presented. The time-space 
relation of the drift chambers and the procedure for the time-of-flight calibration is de- 
scribed. Also, the procedure used for the monitoring of the relative beam-target setting is 
discussed. 
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In chapter [5] the method of the pp — > pp-q' reaction identification will be depicted and the 
method of identifying measured and also unobserved particles will be given. 
Further on, the procedure of the luminosity (L) determination will be presented in chapter 
El and in the consecutive chapters it will be shown how the absolute value and the spread of 
the beam momentum was extracted, and the method for the determination of the position 
of the drift chambers relative to the COSY-11 dipole will be explained. 
Chapter [9] comprises the evaluation procedure of the differential distributions. First, the 
kinematical fit procedure and then the background subtraction due to the multi-pion pro- 
duction will be discussed. 

The final results concerning the total cross section and differential distributions are pre- 
sented in chapter [TOl The acceptance corrections are discussed and the achieved experi- 
mental results are compared to theoretical predictions. 
The conclusions are presented in chapter [TTl 

In the appendices at the end of the dissertation some issues discussed in the thesis are 
explained in more detail. In the first one the structure of the pseudoscalar meson nonet, 
meson masses and quark structure are presented. The second one is devoted to the descrip- 
tion of the parameterization of the on-shell proton-proton interaction. General remarks 
about the combined analysis of the r]' meson formalism in photo- and hadro-production 
are presented in the third addendum, and the linear energy dependence of the production 
amplitude is discussed in the last one. 
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2. Low energy interaction within a proton-proton-meson 
system 

The interaction of hadrons is the reflection of the strong force between the quarks, and 
provides information about the hadron structure and the strong interaction itself [1]. In 
the framework of the optical model, the interaction between hadrons can be expressed in 
terms of phase-shifts, which in the zero energy limit are described by the scattering length 
and effective range parameters [1]. These variables are quite well established for the (low- 
energy) nucleon-nucleon interaction [18,19], but they are poorly known for the nucleon- 
meson or meson-meson interactions. The estimated real part of the scattering length of 
the 77-proton potential, depending on the method of the analysis and studied region, is 3 
to 10 times [1] larger than for the vr'^-proton scattering (Op^o = 0.13 fm) [20,21], while for 
the rj meson only an upper limit is known of 0.8 fm [22]. 

The interaction of mesons (e.g. pseudoscalar mesono ^r, r], r]') with nucleons could be 
deduced from the experiments realized by means of meson beams, but such experiments are 
not feasible in case of the flavour neutral mesons due to their short lifetime [1,5]. However, 
the study of their interaction with hadrons is certainly accessible via their influence on 
the cross section of reactions like A^A^ — > NN Meson in which they are produced [1]. In 
such a case, the interaction within the final meson-nucleon system will modify the shape 
of the excitation function and of the differential distributions of invariant masses of the 
nucleon-nucleon-meson systems. 

2.1 Excitation functions for pp ppi] and pp pprj' reactions 

Near the kinematical threshold measurements of nucleon-nucleon collisions allow to study 
the particle production with a dominant contribution from one partial wave only. In this 
energy range, the dependence of the total cross section as a function of the centre-of- 
mass excess energy is predominantly determined by the available phase space and the 
interaction between the exit particles. The excitation functions for the pp — > pprj' [23-28] 
and pp — > ppr] [28-33] reactions are presented in figure 12.11 Comparing the data to the 
arbitrarily normalized phase-space integral reveals that proton-proton FSI enhances the 
total cross section by more than one order of magnitude for low energies. In case of the 
rj' meson production one recognizes that the data are described well assuming that the 
on-shell proton-proton amplitude exclusively determines the phase-space population. 



^The features of the pseudoscalar meson nonet are described in appendix [AI 
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1 10 10 

Q [ MeV ] 

Figure 2.1: The pp pprj and pp ppr/' excitation functions [23-33]. The dashed lines indicate 
arbitrarily normalized functions obtained under the assumption of the homogeneous phase space 
occupation. Solid lines correspond to calculations of the phase space weighted by the proton-proton 
on-shell scattering amplitude [22]. 

This indicates that the proton-77' interaction is too small to manifest itself in the excitation 
function within the presently achieved statistical uncertainty. However, for the r] meson 
production the enhancement is by about a factor of two larger than in case of the r]' meson 
and cannot be described by the pp-FSl only. 

2.2 Comparison of 77, 77' and 7r° meson interaction with protons 

The strength of the interaction deduced from the comparison of the data and the lines in 
figure 12.11 depends on the model of the proton-proton interaction used in the calculations 
for the ppr] and ppr]' systems [22]. Therefore, in order to estimate a relative strength be- 
tween the pr] and pr]' interactions in a model independent way one can compare the shape 
of the excitation function of the pp — > ppr] and pp — > ppr]' reactions. Moreover, one can gain 
some quantitative information about these interactions by a comparison of these shapes to 
the ppn^ system [1], since the vr'^-proton scattering length is well known and amounts to 
0.13 fm [20,21]. 
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Figure 2.2: The ratios |Af^|/|MJ| (left) and \M^ |/|Mo I ("ght) extracted from data calculat- 
ing the pp-FSl according to the formulas from reference [22], and neglecting the proton-meson 
interaction [1]. The ratio is shown as a function of the phase space volume [22]. 



For that purpose, one can compare only the dependence of the production amplitudes 
I Mo I derived from the data taking into account the pp-FSl only. 

The dependence of |M^| and | as a function of the phase space volumes normalized 
to \Mq \ are presented in figure [221 The \Mq\ for the ry, rj' and mesons were extracted 
from data, disregarding any proton-meson interaction. When the neglected 77(7?') -proton 
interaction would have been the same as the one for proton- tt", the points in the plots 
should have been consitent with unity as can be seen for the pp — > pprj' reaction, when 
really the interaction shows its weakness, independently of the prescription used for the 
proton-proton final state interaction [22]. In the case of the r/' meson production its weak 
interaction with nucleons at the low-energy range is expected due to the lack of any bary- 
onic resonances which could decay into a Nrj' system [1,8|§. 



Statistical uncertainties allowed to get only a very conservative upper limit for the real 
part of the scattering length of the proton-ry' potential resulting in: 
\Re apr,'\ < 0.8 fm [1,26]. 

Thus, independent of the model used for the prescription of the pp-FSl, from a comparison 
of the energy dependence of the production amplitudes for the pp pprj, pp pprj' and 
pp — > ppvr^ reactions, it was concluded that the interaction within the proton-??' system is 
much weaker than the interaction between the proton and the r] meson. 
Another possibility of learning about interactions within nucleon-nucleon-meson systems 

^Recent calculations of K. Nakayama, Y. Oh and H. Haberzettl predict also resonance state contribu- 
tions for the photo- and hadro-productions [34]. 
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is given by the differential distributions of the invariant masses. This is why the present 
analysis of the pprj' system has been performed in a similar way as it has been done earlier [9] 
for the ppr] system. The determined pp and p-meson invariant mass distributions will be 
used for a comparative study of the interaction within the proton-meson system. 
In the next section the definitions of the studied observables which will be used in the 
further analysis are presented. 

2.3 Definitions of observables 

To describe the studied three particle [ppr]') system one needs only five independent vari- 
ables in the centre-of-mass system. In this frame, due to energy and momentum conserva- 
tion, momentum vectors of protons and rj' lie in one plane, called reaction plane. In that 
plane (shown schematically in figure 12. 3p the relative momenta of particles are described 
by only two variables. These quantities may be chosen as square of the proton-proton in- 
variant mass Spp and square of the proton- ry' invariant mass Sj,^' . Invariant masses depend 
on the relative velocity of the particles and are therefore well suited for a description of 
the interactions between these particles. Besides the relative movement of particles on the 
reaction plane three other variables have to be defined for fixing the orientation of the 
reaction plane in the coordinate system. 



Figure 2.3: Schematic definitions of the centre-of-mass kinematical variables used for the descrip- 
tion of the ppr]' system. In the centre-of-mass frame the momentum vectors of the three outgoing 
particles are located within the reaction plane. In this plane the relative motion of the ejectiles is 
fixed by the square of the invariant masses Spp and Sprj' . Three remaining variables the 0*/ , 6**, 
and ^ angles are used to define the orientation of the emission plane in space. 

In this thesis, by analogy to the evaluation of the ppr] system [9], the azimuthal and polar 
angles of the rj' meson momentum vector relative to the beam direction, denoted as (/>*, 
and 6*1 are used, respectively, and the V' angle describing the rotation of the reaction plane 
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around the direction of the momentum vector of the r/' meson. 

The interaction between final state particles does not alter the orientation of the reaction 
plane [9]. Therefore, it will manifest itself only in the distribution of the invariant masses 
Spp or Sprj', or generally in the population of the Dalitz plot (spp vs. Sp^'). 
In the case of non-interacting particles in the final state these distributions should cor- 
respond to a homogeneously populated phase space. Therefore, their interaction should 
show up as a deviation from these expectation. 

2.4 pp and p — meson invariant mass distributions 

Only two invariant masses of three subsystems are independent and therefore the whole 
accessible information about the final state interaction can be shown in the Dalitz plot. 
One can also use the projection of the phase-space distribution onto the invariant masses 
of proton-proton or proton-meson subsystems [9]. 

The qualitative phenomenological analysis of the determined differential invariant proton- 
proton and proton-r/ mass distributions revealed an enhancement of the population density 
at the kinematical region corresponding to a small proton-r/ momentum. The proton-proton 
and proton-r/ invariant mass distributions determined for the pp pprj reaction at an ex- 



of the figure depict the results of calculations where only the on-shell amplitude of the 
proton-proton interaction has been taken into account. 

In those calculations the enhancement factor has been estimated as the square of 
the on-shell proton-proton scattering amplitude derived using the modified Cini-Fubini- 
Stanghellini formula including the Wong-Noyes Coulomb corrections [22]. 
One can easily see that the mentioned effect is too large to be described by the on-shell 
inclusion of the proton-proton FSI. 

In fact a better description is achieved when contributions from higher partial waves 
or off-shell effects of the proton-proton potential are taken into account. These calcula- 
tions compared to the experimentally determined differential proton-proton invariant mass 
distribution are presented in figure 12.51 In the left panel of this figure the experimentally 
determined differential cross section as a function of the squared invariant proton-proton 
mass is compared to the calculations of V. Baru and collaborators [39] under the assump- 
tion of a ^Pq -^^Sqs transition according to the models described in [39], depicted as the 
solid line. 

Dashed and dotted lines on the left panel of figure 12.51 represent the calculations of 
K. Nakayama and his group [12]. The authors claim that the contribution of the S-wave 
alone is unable to explain the observed enhancement in the squared proton-proton invariant 

similar enhancement has been observed for several measurements performed at excess energies of 
4.5 MeV [9], 10 MeV [35,36], 15 and 41 MeV [37], and 72 MeV [38]. 



cess energy of 15.5 MeV are presented 




The dashed lines in both panels 
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Figure 2.4: Distributions of the square of the proton-proton (spp) (left) and proton-77 (sj,,,) (right) 
invariant masses determined experimentally for the pp — > pp-q reaction (closed squares) . The inte- 
grals of the phase space weighted by the square of the proton-proton on-shell scattering amplitude 
(dotted lines)-FSIpp have been normalized arbitrarily at small values of Spp. The expectations 
under the assumption of a homogeneously populated phase space are shown as solid lines. 
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Figure 2.5: (Left panel) Distribution of the square of the proton-proton (spp) invariant mass 
determined experimentally for the pp ppr/ reaction. Solid and dashed lines correspond to 
calculations under the assumption of a ^Pq -^^Sqs transition according to the models described 
in [39] and [12], respectively. The dotted line shows the result of calculations with inclusion of the 
^5*0 ^^Pqs contribution as suggested in [12]. (Right panel) The same data as in the left panel 
but compared with three-body calculations [10,11]. The solid line was determined with a rigorous 
three-body approach. The dashed line depicts the situation if only pairwise interactions {pp + prj) 
are allowed. 
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mass distribution. Seeking for the better description they postulate that the shape of the 
enhancement can be reproduced by folding the relative momentum of the proton-proton 
subsystem with the available phase space [12] suggesting that the enhancement could be 
the consequence of the pp P-wave in the final state. 

Calculations assuming a ^Pq — >^ Sqs transition correspond to the dashed curve and the 
result of calculations with the inclusion of the ^Sq -^^ Pqs contribution is depicted by 
the dotted line. Although the dotted line corresponding to calculations based on the 
stronger P-wave contribution is in quite good agreement to the experimental determined 
differential distribution of the proton-proton invariant mass, it underestimates the to- 
tal cross section data taken for the pp — > pprj reaction near the kinematical threshold 
(Q lower than 30 MeV) [1]. 

On the other hand, the discussed effect can in principle be assigned to changes of the 
production amplitude, since in calculations by V. Baru et al. [39], and by K. Nakayama 
and collaborators [12] the production amplitude was nearly constant. 
An analysis guided by the assumption of a linear energy dependence of the production am- 
plitude was performed by A. Deloff [13]. The squared invariant mass distributions Spp and 
Sprj determined for the pp —> ppr] reaction measured at an excess energy of Q = 15.5 MeV 
are compared to calculations performed by A. Deloff in figure [2?6l 




Spp[GevV] Sp^[GeVV] 

Figure 2.6: Distributions of the square of the proton-proton (spp) (left) and proton-?/ (sp^) (right) 
invariant masses determined experimentally for the pp pprj reaction (closed squares). The 
experimental data are compared to calculations performed assuming the linear energy dependence 
of the production amplitude as proposed by A. Deloff [13] depicted by solid lines. 

The squared invariant mass distributions could be quite well reproduced by lifting the 
standard on-shell approximation in the enhancement factor and allowing for a linear energy 
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dependence in the leading ^Pq -^^Sqs partial wave amplitude [13]. Those calculations are 
in contradiction to the suggestion of Nakayama [12], giving evidence that higher partial 
waves play only a marginal role. 

At this point, the observed enhancement could be explained by three different hypothe- 
ses: 

i) a significant role of proton-r/ interaction in the final state, 

ii) an admixture of higher partial waves or 

iii) an energy dependence of the production amplitude. 

Based on the pp pprj data only, it is not possible to verify any of those models. These 
contingencies motivated the work presented in this thesis which is an analysis of a high 
statistics pp — > pp-q' reaction measurement in order to determine the distribution of events 
over the phase space for an excess energy of Q = 15.5 MeV the same one as selected before 
for the pp pprj reaction. The comparison of the differential distributions for the proton- 
proton and proton-meson invariant masses in the r/ and 77' production could help to judge 
between postulated explanations of the observed effect and may allow for a quantitative 
estimation of the proton-ry and proton-ry' interaction. 

The experimental facility, the method of the analysis, and achieved results for the pp — > ppr]' 
reaction will be presented in the following chapters. 



3. Experimental facility 



The measurement of the pp — > pprj' reaction was conducted using the cooler synchrotron 
COSY and the COSY-11 detector setup. Both facihties will be described in this chapter. 



3.1 Cooler Synchrotron COSY 

The COoler SYnchrotron (COSY) [17] is located at the Institute of Nuclear Physics of the 
Research Centre Jiilich in Germany. The facility was designed to accelerate polarized and 
unpolarized proton and deuteron beams in the momentum range from 0.3 GeV/c up to 
3.7 GeV/c. The sketch of the whole accelerator complex is presented in figure ISTTl The 
total length of the synchrotron ring is 184 meter. There are two straight 40 meter sections, 
and two bending sections with 24 dipole magnets. 

The experimental installations at the synchrotron can be classified as two groups; 

a) the detectors installed inside the COSY ring: WASA [40,41], COSY-11 [14-16], PISA [42], 
EDDA [43], COSY-13 [44], and ANKE [45] and 

b) outside of the COSY ring at external beam lines: COSY-TOF [46], JESSICA [47], 
NESSI [48], GEM [49], MOMO [50], and HIRES [51]. 

Some of those experiments are already completed and no longer in operation (labelled in 
black) and the others are still in operation (labelled in green) in figure ISTTl 

The COSY synchrotron is equipped with electron and stochastic cooling devices which 
are used to decrease the momentum and spatial spread of the beam [52]. 
In the case of electron cooling, the electrons with velocities equal to the nominal proton 
beam velocity are injected at a straight section of the synchrotron. This operation causes 
that faster protons are decelerated and slower ones are accelerated. 

Stochastic cooling uses an electromagnetic device, the so called pick-up unit, which mea- 
sures the beam deviation from the nominal position at one point of the accelerator and 
corrects it by transmitting a correction signal through the shortest way to the kicker 
unit at the other side of the beam pipe. It causes not only a shift to the nominal beam 
orbit, but also decreases the spread of the transversal and longitudinal momentum com- 
ponents [53, 54]. 
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Figure 3.1: Schematic view of the floor plan of the COSY synchrotron. Marked in violet are 
internal-beam [14,40,42-45] and external-beam [46-51] detector setups. Facilities COSY-TOF, 
WASA and ANKE labelled in green are presently still in operation. The positions of electron and 
stochastic cooling devices are indicated and the stochastic cooling signal lines are presented by the 
light-blue dashed lines. 
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3.2 COSY-11 detector setup 

The measurement of the pp — > pprj' reaction is based on the registration of the two outgoing 
protons and reconstruction of their momenta. The ij' meson is identified using the missing 
mass technique. 




Figure 3.2: Schematic view of the COSY-11 detection facility [14]. Note, that only those detectors 
which were used during the measurement of the pp ppX reaction are presented. Protons 
originating from the pp — > ppX reaction are bent in the dipole magnetic field, and leave the 
vacuum chamber through the exit window. Afterwards they are detected in the two drift chambers 
Dl and D2, in the scintillator hodoscopes SI and S2, and in the scintillator wall S3. The scintillation 
detector S4 and the silicon pad detector Si are used in coincidence with the Dl, D2 and SI detectors 
for the registration of the elastically scattered protons. 
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The COSY-11 facility is one of the internal detector setups installed inside the COSY 
synchrotron tunnel at a bending section of the ring. It is mounted next to one of the 
dipole magnets, and benefits from the dipole magnetic field which is used for the particle 
separation from the beam. 

A schematic view of the COSY-11 apparatus is presented in figure 13.21 The figure 
illustrates also schematically the tracks of protons outgoing from the pp ppX reaction. 
Two outgoing protons possessing smaller momenta than the beam momentum, are bent in 
the dipole magnetic field towards the detector system. They leave the vacuum chamber 
throughout the exit window made out of a 30 /im layer of aluminum and 300 /-tm of a 
carbon fiber carrier material with an averaged density of 2.1 g/cm^ [14] and are detected 
using the drift chambers Dl and D2, the scintillator hodoscopes SI and S2, and the scin- 
tillation wall S20. 

The targeJnl used during the experiment, was realized as a beam of H2 molecules 
grouped inside clusters of up to about 10^ atoms. The average density of the target 
was around 5 • 10^^ atoms/cm^ [63]. It was installed in front of the dipole magnet as it can 
be seen schematically in figure 13.21 

The drift chambers Dl and D2 were used for the determination of the particles trajec- 
tories. Those two planar drift chamber stacks are spaced by 70 cm [14,64]. Their active 
area is 1680 mm wide and 433 mm high. Drift chamber Dl (standing closer to the bending 
magnet) consists of six detection planes. The first two with vertical wires, two with wires 
inclined by +31° and two inclined by —31°. The D2 drift chamber is built in the same 
scheme, but it is extended by two additional planes with vertical wires. 

The wires in adjacent planes of each pair are shifted by half of the cell width to resolve 
the left-right position ambiguity with respect to the sense wire. The chosen configuration 
of the detection planes allows to perform the measurement of the horizontal and vertical 
coordinates and enables a unique multi-hit event identification [64]. 

A charged particle crossing the drift chambers produces gas ionization inside the drift 
cells, filled with a gas mixture of one to one argon and ethane at atmospheric pressure. 
The electron drift time to the sense wire is a measure of the distance between the passing 
particle track and the sense wire (see section [4?T]) . In the case of particle tracks oriented 
perpendicular to the detection planes, the maximum drift time corresponding to the max- 
imum drift path of 20 mm equals to 400 ns. 



^In the measurement of the meson production in the quasi-free pd pnpspectatorX reaction dedicated 
neutron [55-58] and spectator [59-61] detectors were installed in addition. They were however not used 
for the measurement described in this thesis. 

^The H2 cluster target specifications are described in the references ]62, 63]. The dimensions of the 
cluster target used during the measurement are described in details in chapter [7] 
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Determined particle trajectories in the data analysis are traced through the magnetic 
field of the dipole back to the target. Therefore, it is possible to reconstruct the momentum 
vectors of outgoing particles at the reaction point. The reconstruction of the momentum 
vectors of the registered particles combined with the information about the time-of-flight 
between the SI and the S3 detectors allows for the calculation of the particle mass and by 
this the particle type identification. 

The SI scintillating hodoscope is built out of sixteen identical, vertically arranged 
modules, read out from both sides (top and bottom) by photomultipliers. The modules 
with 45 X 10 X 0.4 cm^ dimensions are arranged with small vertical overlap (1 mm [14,65]) 
in order to avoid "not covered" space in the geometrical acceptance. The SI detector is 
used as the "start" for the time-of-flight measurement. 

The S2 scintillating hodoscope, similar as SI, consists of sixteen scintillation modules 
with the dimension of 45 x 1.3 x 0.2 cm'^ [14]. 

The S3 scintillating detector delivers the "stop" information for the time-of-flight mea- 
surement. It is built from one non-segmented scintillating wall with the dimension of 
220 X 100 X 5 cm^. It is viewed by a matrix of 217 photomultipliers [14,66,67], occupying 
the edges of equilateral triangles with the sides of 11.5 cm. 

The S4 scintillation counter together with the silicon pad detector (depicted in figure [32] 
as Si) are used for registration of the recoil protons from the proton-proton elastic scattering 
[14]. The silicon pad detector [14] consists of 144 pads with dimensions of 22.0 x 4.5 x 
0.28 mm^. Each pad is read out separately. 

3.3 Trigger logic 

In the experiment two independent trigger branches were used, in order to detect the 
pp — > pprj' and pp pp reactions. 

The main trigger used for the pp pprj' reaction was based on the following conditions: 

Tpp^pprf' = {Sl^>2 V Sl^=i^high V S2^>2 V 32^=1 high) A 5'3^>2, (3.1) 

where fj, denotes the multiplicity of segments in the SI and S2 scintillation hodoscopes, and 
the number of photomultipliers which have fired in the S3 detector. The subscript high 
stands for a high amplitude signal in the SI and S2 detectors which was implemented for 
triggering events when two particles cross the same segment [67]. The hardware threshold 
for the high amplitude was set high enough to reduce the number of single particle events 
considerably, and low enough to accept most events (almost 100 %) with two protons 
passing through one segment of SI [67]. 

The trigger used for the selection of elastic scattering reactions was based on the coincidence 
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between signals from the SI and S4 detectors: 

Tpp^pp = 5'l^=i A 34, (3.2) 

where the SI hodoscope was used for the registration of forward scattered protons and the 
S4 scintillation detector was registering recoil protons. 

The detectors were positioned to cover a large part of the kinematics of the pp pp 
elastic scattering. Due to the high rate of the pp — > pp reaction only every 128'th event 
was registered for the further analysis. The number of the elastic scattering events was 
later used for the luminosity determination (see chapter [6|) . 



4. Calibration of the detector setup 



In this chapter the method used to cahbrate the COSY- 11 detectors and their relative 
settings will be presented. In particular, the time-space relation of the drift chambers and 
the procedure of time-of-flight calibration will be described. In addition, the procedure of 
monitoring the relative beam-target setting will be discussed. 



4.1 Space-time relation for drift chambers 

The drift chambers Dl and D2 consist of 6 and 8 planes of wires, respectively. They provide 
the information about the drift time of electrons (to the sense wires) produced along the 
trajectory of charged particles passing through the chambers. In order to reconstruct those 
trajectories one needs to establish a relation between drift time and distance between the 
particle track and the sense wire (Fig. 14.11 (left)). 
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Figure 4.1: (Left panel) Time-space calibration - relation between the distance from the sense 
wire and the drift time. (Right panel) Corrections of the time-space relation, where the middle 
histogram represents the mean value of AX, the upper and lower histograms visualize one standard 
deviation (tr) of the AX distribution. See the text for the explanation of AX. 



Due to the sensitivity of the drift velocity to the atmospheric pressure, humidity and 
gas mixture changes [68], the data used for the calibration procedure were divided into 
time intervals of about 3-8 hours, with a similar number of collected events. The cal- 
ibration function was determined for each interval applying the procedure of iterative 
improvements [67]. Starting with an approximated function of the space-time relation, the 
distance X from the particle trajectory to the sense wire has been calculated. Then, a 
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straight line to the obtained points was fitted and further on, assuming that it corresponds 
to the real particle track, the deviation AX , between the measured distance of the particle 
track from the sense wire and the one from the fitting procedure was calculated. Next, 
having a certain amount of data, one could determine a mean value of AX as a function 
of the drift time (presented in figure |4?T] (right)). The /S.Xmean was subsequently used 
for a correction of the time-space relation. Next, the improved function has been used for 
the track reconstruction similarly as in the first step. The whole procedure was repeated 
until the corrections were smaller than the statistical uncertainty of the /S.Xmean- 
The averaged spatial resolution of the drift chambers achieved in the experiment discussed 
in this thesis amounted to 250 /im (rms). 

4.2 Time-of-flight calibration of scintillator detectors 

The scintillator detectors SI and S3 "start" and "stop", respectively, are used for the time- 
of-flight measurement. SI consists of 16 scintillator plates with photomultiplier readout 
from both sides, and S3 is a scintillator wall read out by a 217 photomultipliers matrix. 
In order to obtain the proper information about the time-of-flight between both detectors, 
one needs to determine time "offsets" for all photomultipliers i.e. the relative differences 
in transition time of the signal from the photomultiplier to the TDC unit. 
Let us denote tsi as the real time when a particle crosses the SI detector and ^53 when it 
crosses S3. Then, the time-of-flight can be calculated as follows: tofsi~s3 = tss — tsi- 
The measured TDC values for a single photomultiplier in SI and S3 detectors read: 

TDCsi = tsi +ty+ Cai,{PM) + tl}j^^,{PM) - ttrigger, (4.1) 
TDC S3 = ts3 + tpos + fllikiPM) + t'!]j^^t{PM) - ttrigger- (4.2) 

In both equations the time stamp ttrigger (denoting the time of the trigger signal) is iden- 
tical. The index y corresponds to the distance between hit position and the edge of the 
scintillator close to the given photomultiplier in the SI detector and pos stands for the 
distance between the hit position and the photomultiplier in the S3 detector. The ab- 
breviation tyjaik denotes the corrections for the time — walk effect, i.e. the signal time 
dependence on the signal amplitude [69]. Any dependence of ty is cancelled by taking the 
average between the times measured by the upper and lower photomultipliers [67] and this 
can be calculated from the known trajectories. Thus, the only unknown variables are the 
time offsets io//set for both detectors. For a first approximation, the time difference in 
the SI detector can be achieved by taking into account signals from the particles cross- 
ing the overlapping parts of the modules. Next, for the S3 detector the time offset can 
be extracted from the comparison between tofsi-S3 = ts3 — tsi and the time-of-flight 
calculated from the particle momentum reconstructed via curvature in the magnetic field 
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{to free)- Then, iteratively, using the obtained S3 offset one can determine the time off- 
sets for the SI detector. After two iterations, the time offsets for both detectors can be 
estabhshed. As an example, the distribution At{PMj^) determined as 

At(PMfJ) = to/5i_s3(PM/^) - tofrec 

for a group of photomultiphers of the S3 detector are presented in figure 14.21 
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Figure 4.2: Distribution of the time difference between the time tofrec calculated from the re- 
constructed particle momentum and the measured time tofsi-s3 between the SI detector and a 
particular photomultipHer in the S3 detector, as a function of the photomultipHer ID in the S3 
detector. As achieved after second iteration, the figure depicts only a fraction of PM's of the S3 
detector (75-112). 

The time offsets for the photomultiphers in the S3 detector are obtained on the basis of the 
time differences between tofsi-S3 and tofrec presented in figure 14. 2[ They were adjusted 
such that this difference is equal to zero. 

4.3 Monitoring of relative beam-target settings 

Possible changes of the position where the beam crosses the target could have significantly 
infiuenced the momentum reconstruction and as a consequence could worsen the determi- 
nation of the mass of the undetected particle. Therefore, it is important to monitor the 
position of the beam and target overlap. The center of the beam-target overlap can be 
determined from the momentum distribution of the elastically scattered protons [70]. 
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Figure 4.3: The distribution of the perpendicular pj_ versus parallel P|| momentum components for 
pp pp elastic scattering events at a beam momentum of 3.260 GeV/c. The solid Hne corresponds 
to the kinematical ellipse. Changes of the event density along the kinematical ellipse reflect the 
angular dependence of the cross section for the pp — *■ pp reaction. 

The mean value of the distance between the expected kinematical ellipse and the exper- 
imental points (shown in figure 14. Sp may be used as a measure for the deviation of the 
center of the interaction region from its nominal position (Acenter). A pictorial definition 
of Acenter is presented in figure 14.41 and the beam-target geometrical conditions are de- 
picted in figure 14.51 By assuming a wrong interaction center the reconstruction results in 
a wrong momentum determination and the pp — > pp events are not centered around the 
expected kinematical ellipse. 
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Figure 4.4: Pictorial deflnition of the deviation of the center of the interaction region from its 
nominal position (Acenter). 
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target 






Figure 4.5: Schematic view of the relative target and beam settings. Left panel depicts the view 
from above, right presents a side view, ux and cry denote the horizontal and vertical standard 
deviations of the assumed Gaussian distributions of the proton beam density, respectively. Ax^j 
denotes the distance between the centers of the proton target and beam. The figure is adapted 
from [70]. 

In the left panel of figure 14.61 the mean distance of experimental pp pp events from the 
expected kinematical ellipse is shown as a function of Acenter assumed in the analysis. As 
can be seen, the center of the interaction region differs by 0.45 cm from the nominal one. 
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Figure 4.6: (Left panel) The distance between the expected elHpse and the center of experimental 
distribution on the (p_L,p||) plot (figure [431) versus the deviation of the center of the interaction 
region from its nominal position (Acenter) . (Right panel) The deviation of the distance from the 
kinematical ellipse as a function of the time of the measurement. The mean value of the distance 
from the ellipse has been plotted for 13 hours intervals. In the analysis the value of /^center was 
set to 0.45 cm. 



In the right panel of figure [461 it is demonstrated that the beam-target conditions were 
stable during the course of the experiment. Fluctuations seen in the figure are within the 
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statistical error in the determination of the mean value of the distance from the ellipse. 
The variations are at a level of 10~^ and, as can be inferred from the plot presented in the 
left panel of figure [46l correspond to shifts of the interaction center by less than 0.01 mm. 
Thus, the variations of the center of the interaction region can be safely neglected in the 
further analysis. 



5. Identification of the pp pprj' reaction 

In the following chapter the method of identifying the pp pprj' reaction will be described. 
5.1 Identification of protons 

The measurement was based on the registration of two outgoing protons originating from 
the pp — > ppX reaction. 
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Figure 5.1: Distribution of the squared invariant mass of the registered particles. (Note the 
logarithmic scale on vertical axis.) Signals from measured pions, protons, and deuterons are easily 
distinguished. 



After choosing only two-track events, the protons were identified by the determination of 
their rest masses. The particle mass was calculated according to the following formula: 



m 



(3' ■ 



(5.1) 



where p and /? are denoting the momenta and velocities of particles, respectively, which 
were determined in an independent way {p from the curvature of the trajectory in the 
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dipole and (5 from the time of flight between SI and S3). The distribution of the squared 
masses of the particles is shown in figure ISTTl Clearly visible are signals from pions, protons 
and deuterons. 

For the further analysis, particles with reconstructed masses in the range from 0.2 to 
1.5 GeV^/c^ were assumed to be protons. 

5.2 Identification of the ri' meson 

In the present experiment the decay products of the rj' meson were not measured, therefore 
it was impossible to identify its production on an event-by-event basis. Even in experiments 
detecting all decay products an unambiguous identification of a single "rj' production event" 
is not possible, but the background would be much smaller. 
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Figure 5.2: Missing mass spectrum for the pp — > ppX reaction measured at a beam momentum 
of Pb = 3.260 GeV/c. The peak originating from the pp — > pp-q' reaction is clearly seen on top 
of a multi-pion production background. The dashed line corresponds to a fit with a second order 
polynomial to the data outside the signal from the rj' meson. 

The number of pp pprj' events was determined using the missing mass technique. This 
method is based on the knowledge of the protons four-momenta before and after the re- 
action. Denoting: P;, = {Eb,Pb), = {Et,0), Pi = and P2 = {E2,P2) as the 



Identification of the pp pprj' reaction 



37 



four-momenta of the proton beam, proton target, and first and second outgoing proton, 
respectively, one can use the following formula, in the case of the pp — > ppX reaction, to 
calculate the mass mx of the unregistered particle: 

= {Eb + Et-Ei- E2f - {pb +Pt-Pi- nf. (5.2) 

In figure 15.21 the missing mass spectrum determined experimentally for the pp — > ppX 
reaction for the whole data sample is presented. The spectrum includes a broad distribu- 
tion from multi-meson production and the well defined peak originating from the r]' meson 
production. 

The smooth behaviour of the experimental multi-pion production background, which could 
be verified by Monte Carlo simulations studies (see section [9?2]) . allows for a simple poly- 
nomial fit. The knowledge of the smooth behaviour of the cross section [60,71], assures 
that in the range of the signal, the multi-pion background should be flat. 
The dashed line in figure 15.21 corresponds to a second order polynomial fitted to the ex- 
perimental background. Indeed, it can be seen that the fit reproduces the shape of the 
background satisfactory well. 
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Figure 5.3: Result of the background subtraction using the second order polynomial fit presented 
in figure [531 as a dashed line. 
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In figure 15.31 the experimental missing mass spectrum after the background subtraction is 
presented. The background was approximated by a second order polynomial fit as depicted 
in figure E2] by a dashed line. The total number of registered and reconstructed pp — > pprj' 
reactions amounts to about 15000. 

Here, the statistics achieved in the measurement is only illustrated and the possibility of 
the background determination is shown. A detailed discussion of the subtraction of the 
multi-pion production background for differential cross sections will be comprehensively 
described in section [921 



6. Luminosity determination 



In order to determine the absolute values of the differential cross sections, the luminosity 
(L) integrated over the measurement time has to be established. For that purpose, the 
analysis of the pp — > pp reaction, in order to establish the number of elastic scattering 
events was performed. 

A schematic view of the COSY-11 detector setup with superimposed tracks of elastically 
scattered protons is shown in figure 16.11 




Figure 6.1: Schematic view of the COSY-11 apparatus presenting the detectors used for the reg- 
istration of the pp pp elastic scattering. The superimposed lines show an example of trajectories 
from elastically scattered protons with the laboratory angles 9i and 02- One can compare this 
figure to figure [521 The xsi axis indicates the size of detector modules (for further description see 
text). 
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One can evaluate the luminosity (L) according to the formula: 

AiV(^2) _ da* 

where ^^(^2) denotes the known differential cross section [72] and AA^(02) indicates the 
number of elastically scattered protons at a solid angle Ail* around the proton emission 
angle 62 in the centre-of-mass system. In the further analysis, the available range of the 
62 angle (44° to 66°) was divided into 11 bins with a width of 2°. The tracks of elastically 
scattered protons resulting in signals in the SI detector with a coincident signal in S4 from 
the second proton covers the horizontal axis of the SI detector, marked in figure 16.11 as 
xsi, from 40 cm to 75 cm, which corresponds to a 62 angle range from 44° to 66°. 
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Figure 6.2: Projection of the event distribution along the kinematical ellipse for a centre-of-mass 
proton scattering angle of 62 in the range from 48° to 50°, corresponding to the range in the SI 
detector from .T51 = 67.0 cm to xsi = 70.0 cm. 

As an example, the distribution of elastically scattered protons at the O2 angular range 
from 48° to 50° is presented in figure 16.21 This distribution shows the projection of the 
experimental data along the kinematical ellipse. The number of events (reduced by the 
background indicated by the dashed line) is used for the calculation of the luminosity. 
The signal from elastically scattered protons can be clearly separated from the flat multi- 
pion scattering background. The solid angle AfJ* or the certain angular bin, is calculated 
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using the Monte-Carlo method, as follows: 

AO* = ^accepted j^^j^ 2) 

2 A^o 

where A'o stands for the number of proton-proton elastic scattering events in the corre- 
sponding angular range and Naccepted constitutes the number of events in the considered 
bin of the 62 angle, which could be registered and identified. In particular, the analysis in 
the following manner was done. First, A'q = 2 • 10^ events has been generated, calculating 
the response of the COSY-11 detectors, and then those events have been analysed using 
procedures applied for the experimental data evaluation in order to determine the number 
of Naccepted cvcnts for each ^2 angle interval. 
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Figure 6.3: Differential cross section for the proton-proton elastic scattering. The result of this 
thesis (closed circles) measured at a beam momentum of Pb = 3.260 GeV/c was scaled in amplitude 
to the cross section measured by the EDDA collaboration shown by open squares [72]. 

Figure 16.31 indicates the angular distribution of the differential cross section for elastic 
proton-proton scattering obtained in the experiment (closed circles). The amplitude of 
that distribution was fitted to the data of the EDDA experiment including only one free 
parameter being the integrated luminosity (see eq. 16. ip . The extracted integrated lumi- 
nosity for the experiment described in this thesis amounts to L = (5.859 ± 0.055) pb"-*^. 
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The knowledge of the luminosity value will allow for the overall normalization of the 
derived differential cross section as a function of the Spp and Sp^' invariant masses, angular 
distributions and total cross section which will be discussed in chapter [TOl 



7. Determination of the spread and the absolute value of the 
beam momentum 



In order to perform realistic simulations of the studied reactions, in particular to determine 
the acceptance and to calculate the covariance matrix, it is mandatory to know the absolute 
value and the spread of the beam momentum. The discussed measurement of the pp — > pp-q' 
reaction was nominally performed at the same value of excess energy Q as the pp — > pp 
reaction measurement with Q = 15.5 MeV which corresponds to a nominal proton beam 
momentum of Pb = 3.257 GeV/c. 

The precision of the absolute beam momentum adjustment of the COSY synchrotron is 
about 10~^ [52] which in this case corresponds to ~3 MeV/c. 

The beam momentum dependence of the mean value of the missing mass distribution 
presented in figure 15.31 was studied in order to determine the actual value of the beam 
momentum more accurate. 

The beam momentum was calculated using the formula: 

mx = Vs -2mp = (2ml + 2mp^Jp^'+m^^ ^ - 2mp, (7.1) 

where -y/i denotes the total energy in the centre-of-mass frame, Pb stands for the proton 
beam momentum, and nip corresponds to the proton mass. 

The beam momentum of Pb = 3.260 GeV/c has been determined by adjusting the Pb 
such that the mean value of the missing mass peak is equal to the rj' meson mass. The 
determined value of the beam momentum differs by 0.003 GeV/c from the nominal one. 
This deviation is in agreement with results of analogous analysis performed in previous 
measurements [67]. 

The determined value of the excess energy amounts to (Q = 16.39 ± 0.01 ± 0.4) MeV, 
where the errors indicate statistical and systematic uncertainty, respectively. The domi- 
nating systematic uncertainty was established in [9] and the statistic uncertainty of the 
excess energy was determined using the following formula: 

where APb was obtained from the linear relation to the missing mass as described in [67]. 
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Figure 7.1: Spectrum of the beam momentum distribution integrated over the whole measurement 
time. The value of APb = corresponds to a beam momentum of 3.260 GeV/c. The dashed lines 
mark the beam momentum dispersion for the extended 9 millimeter target used in this experiment. 

After the determination of the real value of the absolute beam momentum, now its spread 
will be determined. One can calculate the beam momentum spectrum from the frequency 
spectrum of the COSY beam (Schotky spectrum measured during the experiment) using 
the below formula [73]: 



where / and Pb denote beam frequency and beam momentum, respectively, rjB is a pa- 
rameter which depends on the beam optics, i.e. the electric and magnetic fields in the syn- 
chrotron. During the experiment, the rjs parameter was established to be r]B = 0.12 [74]. 
The spectrum of the beam momentum obtained during the experiment is shown in fig- 
ure [LH 

The dispersion of the beam momentum depends on the magnetic field along the ring. Pro- 
tons on the outer routes have a "longer way" than those on the inner side of the beam, 
and the trajectories are different from the nominal value. When at a certain point the 
particle position deviates by Ax from the nominal (xq) position and possesses the relative 
momentum deviation APb/Pb, one can relate this value by means of the known dispersion 




APb 



Pb 



(7.3) 
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(D) using the following formula: 



where D = f3-Dmod with Dmod = 14 m [74,75] and the particle velocity f3. Applying in the 
upper formulas the value of the beam momentum Pb = 3.260 GeV/c, f3 = 0.961, and 
using a 9 millimeter wide target, one obtains APb = 1-1 MeV/c. This value is depicted 
in figure [7?T] as the area between the dashed lines. Thus, the maximum beam momentum 
spread in the interaction point is APb = ±1-1 MeV/c. 
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8. Fine tunning of the relative dip ole- chamber settings 



In section |4?T] a space-time calibration of the drift chambers was described. In this chapter 
a procedure for the determination of their position relative to the COSY-11 dipole will 
be shown. For this purpose, the missing mass method for the pp — > ppX and the kine- 
matical ellipse for the pp — > pp reaction, as described in the previous chapters, will be used. 

The proper knowledge about the position of the drift chambers is needed in order to 
achieve the optimal resolution for the momentum reconstruction. In the horizontal plane 
the position of the drift chambers stack can be defined by two variables: the horizontal 
shift of the chambers Ax and their inclination Aa. 



Figure 8.1: Schematic view (from top) with the pictorial definitions of the horizontal shift of the 
chambers Aa; and their inclination Aa. These parameters were used in the analysis to determine 
the proper arrangement of the Dl and D2 drift chambers relative to the COSY-11 dipole. 

The parameters Ax and Aa influence the reconstruction of particle momenta and further 
more indirectly influence the missing mass resolution. In order to fix these parameters, 
an analysis of the proton elastic scattering events was performed, since the distribution of 
the perpendicular p± versus parallel p\\ momentum components for the pp — > pp reaction 
is sensitive to changes of the position of the drift chambers. The distribution of the p± 
versus momentum components for the pp — > pp elastic scattering determined at a beam 
momentum of 3.260 GeV/c, presented earlier in chapter [43l is shown in figure [821 In the 
left panel, the situation when the parameters are properly adjusted is presented. In the 
right panel the situation is shown when the inclination Aa was changed by ~ 0.9°. Studies 
with a variation of Aa in the range from —0.9° to 0.9° and of Ax ranging from 0.0 cm to 
1.0 cm were performed, where simultaneously the p± versus spectrum and the changes 
of the width of the rj' peak in the missing mass spectra were controlled. 
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Figure 8.2: Distributions of the perpendicular pj_ versus the parallel p|| momentum components of 
measured protons for pp — > pp elastic scattering events at a beam momentum of 3.260 GeV/c. The 
solid line corresponds to the kinematical ellipse where the elastically scattered events are expected. 
The distribution is shown for Aa = 0.04° (left), and for Aa ~ — 0.9° (right), respectively. 

The lowest width of the missing mass peak (see Fig. I5.2|) and the best fit of the event 
distribution to the kinematical ellipse for the pp — > pp elastic scattering distribution (see 
left panel of figure [8?2]l has been found for the values of A.x = 0.45 cm and Aa = 0.04°. 
A comparison of the left and right histograms in figure 18.21 shows the sensitivity of this 
distribution to the parameter Aa, which could be determined within an accuracy of better 
than 0.01°. 



9. Evaluation of the differential distributions 



9.1 Kinematical fit 

In order to search for small effects like proton-meson interaction on the population density 
of the phase-space, it is of importance to account for any possible changes of the measured 
distributions due to the finite resolution of the detector system, which may alter the shape 
of the spectrum especially close to the kinematical limit. Therefore, one must include the 
experimental resolution in the theoretical calculations and in order to improve the effect 
of the resolution a kinematical fitting of the data was performed. Both procedures require 
the knowledge of the covariance matrix. 



In order to calculate the covariances and variances between each combination of the 
registered proton momentum components, 2 • 10^ reactions of the type pp pprj' have 
been generated including all arrangements and conditions from the experiment. Then, the 
whole simulated sample of events has been analysed with the same procedures as used for 
the experimental data. In this step for each simulated event the pair of the real (gen- 
erated) protons momenta Pi^gen, P2,gen and the pair of momenta reconstructed from the 
simulated response of the detectors: Pi,rec, ^2,rec are accessible. The available kinematical 
information about an event may be expressed in form of the six dimensional momentum 
vector: P = [pix,Piy,Piz,P2x,P2y,P2z] including the components of the reconstructed mo- 
menta of both protons Pi^rec = \pix,Piy,Piz] and i-2,rec = [P2x,P2y,P2z]- The covariance 
between the i*'^ and j*'* components of P was determined as the average of the product of 
deviations between the reconstructed and generated values. The formula used to establish 
the covariance matrix elements, reads as follows: 



CO'u(i,j) ^ ^ '] i.-^i,gen ^i,recon)(.-^j,gen ^j,recon}^ i.^-^) 
n=l 

where PJ^g^n and P^j-econ stand for the generated and reconstructed values of the i*'^ com- 



n=l 

ponent of the vector P from the n*^ event. 

Due to the inherent symmetries of the covariance matrix {cov{i,j) = cov{j,i)) and since 
the measured protons are indistinguishable, there are only 12 independent values which 
determine the 6x6 error matrix unambiguously [9]. 

The covariance matrix (in units of MeV^/c^) determined for the pp pprj' reaction looks 
as follows: 
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As it was already pointed out in the previous paragraph, 6 variables (2 times 3 components 
of the momentum vectors of the two protons) have been measured in the experiment. It 
was assumed in the analysis, that an event with the missing mass equal to the mass of the rj' 
meson (within the experimental resolution) corresponds to the pp pprj' reaction. Under 
this assumption only five of the kinematical variables are independent compared to six 
measured variables. Therefore, a kinematical fitting procedure can be applied to improve 
the effect of the limiting resolution. Hence, the protons momenta were varied demanding 
that the missing mass of the unregistered particle is exactly equal to the known mass of 
the rj' meson and there has been chosen the momentum vector which was the closest to 
the vector determined from the experiment. The inverse of the covariance matrix was used 
as a metric for the distance calculation. The kinematical fit clearly improves the effective 
resolution, what can be seen in figure 19. 1[ 
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Figure 9.1: Spectrum of the difference between the simulated and reconstructed proton momen- 
tum. The dashed line denotes the spectrum before kinematical fit and the solid line corresponds 
to the situation after the fitting procedure. 
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After the kinematical fit was performed for each event aside from the experimentally de- 
termined momentum vector: P'^^^ = \pTx ^PTy ^pTz^ -^pTx ■'PTy' ^pTz"]-: 
the kinematically fitted momenta: P-^** = \P]^x -iP^^y ^P]^z ^P2x ^P2y -iPiz] were obtained. 
The Pf'* vectors are more precise and assure that the determined distributions do not 
spread beyond the kinematical boundaries. 

For the further analysis of events corresponding to the r/' production, the distribution of 
the from the kinematical fit procedure was checked. In figure 19.21 the distribution 
as a function of the missing mass determined for the experimental momentum vectors is 
presented. Having each event described by the two vectors P^^p and P^^^ and the 
of the kinematical fit, the variables Spp, Sprj', ip and |cos0*/| were evaluated by using the 
kinematically fitted vectors P-^** if the x^^value of the fit for the event was lower than 
1.5. However, if the was larger than this limitation, P^^'p were used, since those events 
correspond to the background with a missing mass lower or higher than the ?/ meson mass. 
For the multi-pion production background the kinematical fit procedure is not justified. 
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Figure 9.2: The distribution as a function of the missing mass. The dashed hne at = 1-5 
visuahzes the cut used in the analysis for the separation between kinematically fitted events (P-^**) 
and events without the fit (P'^^p). 

The evaluation of the differential distributions requires only events corresponding to the 
pp — > pprj' reaction. In order to select the number of those events from the multi-pion pro- 
duction reactions, a missing mass spectrum for each studied interval of the variables: Spp, 
Sprj', tp and |cos0*/| was determined. In the next section, the procedure of the background 
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subtraction applied in the analysis will be discussed. 
9.2 Background subtraction 

In the following section the method of the multi-pion background subtraction is presented 
as used in the analysis. 

For the background free determination of the differential Spp, Sp^/, ip and |cos0*, | distribu- 
tions, first the distribution of the considered variable has been divided into a reasonable 
number of bins, and then for each bin a missing mass spectrum was produced and the 
number of the pp — > ppij' events was calculated for each interval of Spp, Sp^^', ip and 10056**, | 
separately. 

The missing mass distributions include a smooth multi-pion production below and above 
the produced rj' meson mass which allows to estimate the background corresponding to the 
multi-pion creation in the range of the r]' meson signal using a polynomial fitting function, 
as it was already shown in section 15.21 In figure 19.31 examples of missing mass spectra for 
two intervals of the invariant proton-proton mass Spp are presented. 
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Figure 9.3: Examples of experimental missing mass spectra for two chosen intervals of Spp. The 
dotted lines indicate the second order polynomial according to the formula [9 . 21 whereas the dashed 
lines show the result of the fit of two Gaussian distributions described by the formula lOl 

In both panels of figure l9.3l the dotted lines correspond to the polynomial fits described 
by the mathematical function: 

P{mm,a,b,c) = a + h-mm + c-rarn?, (9-2) 

where mm denotes the missing mass and a, h and c are the free parameters varied during 
the fit. Since not only a second order polynomial can be used as a good approximation of 
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the background also a fit by the sum of two Gaussian distributions was performed. The 
dashed hues in both panels of figure [9?3] correspond to the fit with the function G described 
by the following formula: 

G(mm,Ai,Bi,C7i,A2,S2,C2) = Gi {mm, Ai,Bi,Ci) + G2 (mm, ^2, ^2, C2), (9.3) 

where the terms Gi and G2 denote Gaussian distributions depending on the parameters 
Ai, Bi, Gi, A2, B2 and G2 varied freely during the fit and mm stands for the missing 
mass. The fits were performed using the functions of the missing mass (|9.2I and I9.3|l in 
the whole range of missing mass outside of the pp — > ppif signal. As it is seen, the fitted 
functions reproduce the background very well. 

The smooth behaviour of the multi-pion background was verified by Monte Carlo sim- 
ulations. But in order to use these Monte Carlo distributions as a reasonable description 
for the background an extensive time consuming simulation studies would be necessary. 
Therefore, a smooth function adjusted to the regions beside the rj' peak was used instead. 
The sum of two Gaussian distributions was applied in the analysis. It should be stressed 
that the alternative description with the second order polynomial is in rather good agree- 
ment with the result obtained using formula [9.31 

The situation is more complicated for missing mass spectra with the signal close to the 
kinematical limit (see e.g. Fig. 19.41) . In this case the shape of the background on the right 
side of the peak cannot be easily predicted. Such spectra are obtained for kinematical 
regions of higher squared invariant proton-proton masses Spp and relatively low squared 
invariant proton- r/' masses Sp^'. In order to describe the shape of the background in those 
regions, the pp — > pplirr], pp — > pp3iT and pp — > ppAir reaction^] have been simulated and 
the simulated events were analysed in the same way as it was done for the experimental 
data. The result of these simulations is compared to the experimental data in figure 19.41 
The simulations of the different reactions channels were performed with phase space dis- 
tribution including the proton-proton final state interaction [22,60]. 
The simulated missing mass spectra were fitted to the data using the formula: 

5(mm,a,/3,7) = a ■ fpp^pp2nrjimm) + P ■ fpp-,pp3T,{mm) + 'y ■ fpp^ppinimm), (9.4) 

where a, /3 and 7 denote the free parameters varied during the fit procedure. The func- 
tions fpp and fpp->ppATT correspond to the simulated missing mass 
distributions for the pp pp2T:rj, pp — > ppSvr and pp — > pp47r reactions. 



^ Those background reaction channels were chosen as a representation of possible multi-pion production 
background. Since, the pp pp5n and pp pp6n reactions simulations can be neglected and the 
pp ppl-K missing mass of two protons spectrum has similar shape as those for Svr and 47r [71], only 
simulations for listed reactions were performed. 
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In figure 19.41 examples of missing mass spectra for squared invariant proton-proton 
masses Spp € [3.577; 3.580] GeV^/c^ (left panel) and squared invariant proton-r/' masses 
Sp^i € [3.602; 3.605] GeV^/c'^ (right panel) are presented. 
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Figure 9.4: Examples of missing mass spectra determined experimentally for the pp ppX reac- 
tion (solid lines) with superimposed Monte-Carlo simulations (dashed lines) including pp pp2Trr/, 
pp pp^TT and pp — > ppAiT reactions, fitted in amplitude to the experimental data. Missing mass 
spectra for squared invariant proton-proton mass Spp € [3.577; 3.580] GeV^/c^ (left panel) and 
squared invariant proton-77' mass Sprj' G [3.602; 3.605] GeV^/c'' (right panel). The dotted lines in 
both panels correspond to the fit of the formula [9.51 



In both examples, it is clearly seen that the simulations are in a good agreement with the 
experimental background distributions below the r/' peak. Moreover the behaviour of the 
simulated background fits well to the kinematical limit of the missing mass distributions. 
Close to the kinematical boundary, the interaction in the final state is dominated by the 
proton-proton interaction [60] and this interaction influences strongly the shape of the 
observed missing mass spectra. Therefore, it was taken into account in the calculations 
according to the formulas of reference [22] presented in details in appendix iBl 
For the dynamics of the pion production it had been assumed, that pions are produced 
homogeneously over the phase space. As it was described in reference [60], the shape of 
the missing mass spectrum does not change significantly at the edge of the kinematical 
limit if one assumes resonant or direct pion production. 

In order to raise the confidence to the estimation of the background behaviour near the 
kinematical boundary, those distributions were described in an independent way with a 
second order polynomial divided by the Fermi function for the description of the rapid 
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slope at the end of the distributions. To this end the following formula was applied: 

F{mm,a,b,c,d,g) = (a + b ■ mm + c-mm^)/{l + e^™™ " '^^/f), (9.5) 

where a, b, c, d and g are the free parameters varied during the fit procedure. 
The results are presented in figure 19.41 as dotted lines. It is seen that under the rj' peak 
the result of formula [931 agrees well with the background determined from the simulations 
and that both reproduce the shape of the slope quite well. 

In order to perform a further check of the background estimation, examples of the 
missing mass distributions for the regions of the squared invariant masses of proton-proton 
and proton-meson where the rj' is not produced are presented. The missing mass distribu- 
tions are shown in figure 19.51 and represent the regions of low squared invariant masses of 
the proton-meson subsystem (left) and high squared invariant masses of the proton-proton 
subsystem (right). For such values of Spp or Sp^' the production of the r/' meson is not 
possible because Sp^' < (m^ -|- m,,')^ and Spp is too large leaving not enough energy for the 
rj' meson creation. 

The fitted simulations reproduce the background very well. 
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Figure 9.5: Missing mass spectra for low values of Sprj' £ [3.590; 3.593] GeV^/c"* (left) and high 
values of Spp £ [3.586; 3.589] GeV^/c^ (right). The dashed lines correspond to the simulations 
of the pp —* pp2Trr/, pp pp'Sir and pp —* ppAir reactions fitted to the experimental distributions 
(solid lines) using only the amplitudes as free parameters. 
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Figure 9.6: Examples of missing mass spectra determined for the first and fourth bin of the ip 
angle with a width of 0.3 radian. The dashed lines correspond to fits using the function described 
by formula [ 
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Figure 9.7: Examples of missing mass spectra determined for the seventh and tenth bin of the 
ip angle with a width of 0.3 radian. The dashed lines correspond to fits employing the function 
described by formula [931 
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Figure 9.8: Examples of missing mass spectra determined for the first and fourth bin of |cos6'*, | 
with a width corresponding to A|cos0*/| =0.1. The dashed fines correspond to fits according to 
the function described by formula [931 
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Figure 9.9: Examples of missing mass spectra determined for the seventh and tenth bin of \cosd*, \ 
with a width corresponding to A\cos6*,\ — 0.1. The dashed fines correspond to fits using the 
function as described by formula [931 

In case of the angular distributions the background could be described as a polynomial 
or Gaussian in the whole angular range. This is because the signal from the r]' meson is far 
from the kinematical limit for all spectra (see Fig. 19.6119.71 l9.8l and [9?9l) and additionally the 
background varies smoothly on both sides of the peak. In this context it is also important 
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to note that the width of the peak is due to the experimental resolution and since the 
background is smooth on both sides, no structures are being expected below the peak as 
it is verified by the simulation studies including relevant background channels. 
Missing mass distributions for every third bin of ■0 and |cos^*, | are presnted in figures \9M 
19.71 19.81 and 19.91 respectively, using formula [9.31 for fitting to the background. 

9.3 Systematic errors of background estimation 

The main contribution of the systematic uncertainty of the differential cross section deter- 
mination comes from the uncertainty of the estimation of the yield of the r/' events which 
in turn is due to the assumption of the shape of the background. In order to estimate those 
errors, the numbers of background events extracted under two different assumptions were 
compared. For the missing mass spectra with the signal far from the kinematical limit 
the background determined by Gaussian distributions with the background estimated by 
a second order polynomial were compared. Whereas, for the spectra close to the kinemat- 
ical limit a comparison of the background determination by Monte-Carlo simulations with 
a second order polynomial divided by the Fermi distribution was performed. The rela- 
tive difference between determined numbers of events N^^f^ (two Gaussian distributions 
or Monte-Carlo simulations fit, respectively) and N^^j^ (polynomial fit) was used as an 
estimation for the systematic uncertainty of the differential cross sections evaluation. 
The relative systematic errors, expressed in per cent were calculated in the following way: 

_ 

Aback = 'Z' ^ Jb' ■ 100^0- (9-6) 

^^back + ^^back 

The averaged value of the systematic deviations in the various applied fit procedures cal- 
culated according to formula [9.61 are below 3%. 

The systematic uncertainty of the differential cross sections evaluation was calculated using 
the formula: 

'da,,,\ da N^ackU) 



where / stands for both the squared invariant proton-proton mass Spp and the squared 
invariant proton- ry' mass Sp^' . N^ack and A'^^' correspond to the average numbers of back- 
ground and 7]' creation events, respectively. The averaged values of the systematic errors 
for the differential cross sections, calculated according to formula [9.71 are in the order of 
0.32 ;u6/(GeV^/c^) and will be presented in details in the next chapter. 



10. Cross sections 



This chapter is devoted to the derivation of the final results. 



10.1 Acceptance corrections 

The COSY-11 detector setup does not cover the full Att solid angle in the centre-of-mass 
system of the pp — > pp-q' reaction at the proton beam momentum of Pb = 3.260 GeV/c. 
Therefore, in order to study the differential cross sections, one has to perform acceptance 
corrections for the measured distributions. 

In general, the acceptance corrections should be done in the five dimensional phase space 
of e.g. Spp, Spri', (j), ip and cos6*, (as depicted in chapter [2^ . 

But one can safely assume that the dependence on the (j)*, angle, the centre-of-mass az- 
imuthal angle of the -q' meson momentum, is isotropic because of the axial symmetry of 
the initial unpolarized state. In addition, the identical particles in the initial state imply 
that the distribution of cos9*, should be symmetrical around 90°. This limits the variables 
for the differential acceptance to Spp, Sprji , ip and \cos9*,\. 

To determine the distributions of squared invariant masses Spp and Sp^' , the available range 
was divided into 22 bins. Number and width of the bins were chosen as a compromise be- 
tween statistics and the experimental resolution. The width of the bins was chosen for the 
variables Spp and Spr/' to be 0.003 GeV^/c^. 

In the case of angular distributions, the 10056**, | range was divided into 10 bins with a bin 
width of 0.1, and the ip range into 11 bins of 0.3 radian. 

The acceptance correction of the data will be performed iteratively. First, it was assumed 
that the distributions are determined by a homogeneous phase space occupation modulated 
by the pp-FSl. Under this assumption the acceptance was calculated and the differential 
cross section for the variables Spp, Sp^i , tp and |cos^*, | was extracted. Next, using the 
derived distributions the acceptance was calculated and the whole procedure was repeated 
until the differential distributions remain unaltered. 

In order to calculate the differential acceptance it was assumed that the distribution over 
the ^p angle is isotropic like it was experimentally determined for the pp — > ppp, pp — > ppco 
and pp — > ppcj) reactions [76, 77]. 

The geometrical acceptance as a function of Spp and Sprj' is presented in the left and right 
panel of figure llO.il respectively. In figure [T0?2l the acceptance dependence for the |cos6'*, | 
and the ip angular distributions is shown. 
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Figure 10.1: Geometrical acceptance of the COSY-11 detection setup determined for the 
pp ppr/' reaction, simulated for a beam momentum Pb = 3.260 GeV/c, as a function of the 
squared invariant mass of the proton-proton subsystem (left) , and the squared invariant mass of 
the proton-77' subsystem (right). 
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Figure 10.2: The COSY-11 detection acceptance as a function of the polar angle of the 77' meson 
momentum vector |cos6'*/| (left), and for the ip angle describing the rotation of the reaction plane 
around the direction of the ry' meson momentum (right). 



After performing the luminosity determination and having calculated the geometrical 
acceptance, one can evaluate angular distributions. In figure 110.31 the differential distri- 
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bution of the cosine of the 9*, polar angle is presented. The distribution has an isotropic 
character within the statistical errors. 
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Figure 10.3: Distribution of |cos6'*, | of the rj' meson momentum vector in the centre-of-mass 
frame. The shown result was evaluated under the assumption that ^ is isotropic. 

Then, the angular distribution of the tp angle was calculated. Unexpectedly, but similar to 
the situation observed for such a distribution for the pp — > ppr] reaction [9] measured at the 
excess energy of 15.5 MeV, a significant anisotropic behaviour of the ^p angle distribution 
is observed. The ^ distribution determined experimentally is presented in figure 110. 41 
Such a shape cannot be explained by any background behaviour since ^ was extracted for 
each bin separately, as it was mentioned earlier in section 19.21 and in fact the distribution 
shown in figure [TO^ is background- free. Systematic errors in the background subtraction 
are expected to be much smaller and can also not cause such an anisotropy. The missing 
mass spectra used for the background subtraction could satisfactorily well be reproduced 
for each bin in the phase space. For example see figures [9^61 and [9771 in section [921 
Due to the anisotropic behaviour of the ^ the full procedure of acceptance correction 
for the |cos0*, I, Spp and Sp^r distributions was performed once again with the working 
assumption that the distribution ^ is as it was determined from the data. The procedure 
was repeated three times. After all, it was observed that the output distributions are in 
good agreement with the input values. 

It was also checked that the shape of the ^p distribution assumed in the calculations of the 
acceptance does not influence the shape of Spp, Sp^' and |cos0*/| distributions. 
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Figure 10.4: Distribution of the differential cross section as a function of the ip angle de- 
termined experimentally. The superimposed histogram corresponds to a fit by the function 
y = a + b •|sm(-0)|, where a = 7.87 nb/rad and b ~ 1.11 nb/rad. Since, in the analysis the 
advantage of the symmetry ^(V') = H^IV" + was taken, only the range of -0 from to tt is 
presented. 
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Figure 10.5: Distribution of the |cos6'*, | of the 77' meson momentum vector in the centre-of-mass 
frame. The distributions after the second (closed squares) and third (open squares) iterations are 
compared (for details see text). The distributions are nearly identical and therefore the points can 
be hardly distinguished. 
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Figure 10.6: Comparison of the invariant mass Spp distributions after the second (showed as closed 
squares) and third (presented as open squares) iterations assuming different assumption about the 
tp distribution. 
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Figure 10.7: Comparison of the invariant mass Spjji distributions after the second (showed as 
closed squares) and third (presented as open squares) iterations assuming different assumption 
about the ip distribution. 
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As an example figure [T03] presents the \cos6*,\ spectrum determined assuming a ip 
distribution obtained from the second and the third iteration. 

Similarly like in case of the distribution of the cross section as a function of \cos9*,\, the 
differential cross sections as a function of Spp and Sp^' were extracted. First, as it was men- 
tioned, the distributions were determined with the assumption of an isotropic behaviour 
of the angular distributions of polar and azimuthal angle of the r]' meson momentum vec- 
tor in the centre-of-mass frame, and also an isotropic behaviour of the angle distribution. 

Then, again using the experimental parameterization of the distribution of the cross 
section as a function of the 'ip angle in the centre-of-mass frame, the invariant mass Spp and 
Sprj' distributions were determined. Results of the second and third iterations are shown 
in figures 110.61 and 110.71 

Again, it can be observed that the shapes of the distributions after the second and third 
iterations are in agreement. 

10.2 Total and differential cross sections 

In this section, the results of the total cross section and differential cross sections deter- 
mined for the pp — > ppr]' reaction are presented. 

The total cross section for the pp ppr]' reaction, determined at an excess energy of 
Q = 16.4 MeV, is presented in figure [TO^ and amounts to: 139 ±3 nb. In the calculations 
the final distribution of the Spp, Sp^', ip and |cos0*/| were taken into account. 
Though, the shapes of distributions Spp, Sp^' and |cos^*/| are independent of the ip dis- 
tribution, the total cross section depends quite significantly on the shape of ^ [9] which 
is assumed for the acceptance calculations. The obtained value agrees within the errors 
to the previously determined cross section values, however, it is slightly higher than the 
former data. This is due to the fact that in the previous analyses [25-27] the acceptance 
was calculated only approximately disregarding the differential distributions of the ip angle 
which was not established due to lack of statistics. A similar effect of an increase in the 
cross section values after taking into account the differential distribution of ip was obtained 
in the case of the rj meson [9]. 

The systematic error of the total cross section contains a systematic error of the luminosity 
determination which amounts to 3% as it was evaluated in [1] and the uncertainty from 
the ^ distribution, since varying the parameters of the function ^ = a + h ■ \sin{ip)\ 
(within their uncertainties), they change the total cross section by 2.5 nb. The averaged 
systematic error of the background subtraction (like pointed out in section 19. 3|) amounts 
to 5 nb. 
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Figure 10.8: Excitation function for the pp pprj' reaction. The earher measurements [23-28] 
are presented as open squares, the closed black square depicts the total cross section determined 
for the excess energy Q = 16.4 MeV (the result of this thesis). The error bars denote the total 
uncertainty (statistical and systematic). 

Thus, the systematic error of the total cross section adds up to 12 nb. Summarizing, the 
total cross section for the pp pp-q' reaction determined at the excess energy of 16.4 MeV 
is equal to a = (139 it 3 ± 12) nb where the indicated errors are statistical and system- 
atic, respectively. 



In the following tables and figures, the values of the differential cross sections as func- 
tions of the observables |cos0*, |, "0, Spp, and Sp^/ defined earlier in section [2?3l are presented. 
In all tables and figures statistical and systematic errors are included. 
The distribution of the polar angle is listed in table 110.11 and presented in figure 110.91 
The distribution of the ijj angle is listed in table 110.21 and presented in figure 110. lOi 
In tables [TO^ and [TO^ the results for the differential cross sections of the invariant masses 
pp and prj' are presented, respectively. The listed values of both distributions are shown 
in figure 110.111 
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Table 10.1: The differential cross sections as a function of |cos6'*, | for the pp ppij' reaction 
measured at Q = 16.4 MeV. 
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^mmstat 
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0.0024,3,, 



Table 10.2: Differential cross sections as a function of the ■0 angle determined for the pp ppi]' 
reaction measured at Q = 16.4 MeV. Since, ^{ip) = + the full range for ip angle is 

assumed from to tt. 
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Figure 10.9: Distribution of tlie polar angle of the rj' meson emission in the centre-of-mass system. 
The error bars represent statistical and total (statistical + systematic) uncertainties. 
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Figure 10.10: Differential cross section as a function of the angle obtained after the last itera- 
tion. The error bars represent statistical and total (statistical + systematic) uncertainties. 
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Table 10.3: Distribution of the squared invariant mass of the proton-proton system (spp), deter- 
mined for the pp — > ppr]' reaction measured at Q = 16.4 MeV. 
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Table 10.4: Distribution of the squared invariant mass of the proton-77' system (sp,/), determined 
for the pp ppri' reaction measured at Q = 16.4 MeV. 
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Spp[GeVV] Sp^,[GeVV] 

Figure 10.11: Distributions of the squared proton-proton (spp) and proton-77' (sp,,') invariant 
masses, determined experimentally for the pp —* pprj' reaction at the excess energy of Q = 16.4 MeV 
(closed squares). The experimental data are compared to the expectation under the assumption 
of a homogeneously populated phase space which is presented as solid Hues. The integrals of the 
phase space weighted by the proton-proton scattering amplitude - FSIpp are marked as dotted 
histograms. The results of the simulations were normalized to the data in amplitude. 

In figure 110.111 experimentally determined distributions are compared to the expecta- 
tion under the assumption of a homogeneously populated phase space, presented in the 
figure as solid lines. One can see that, for both distributions the data do not agree with 
the calculations. In the same figure the experimental spectra are compared also to the 
integrals of the phase space weighted by the proton-proton scattering amplitude - FSIpp 
which are depicted as dotted histograms. Also the calculations including only the FSIpp 
do not fit to the data satisfactorily. One can see that the comparison depends on the way 
of normalization which will be discussed later. 



10.3 Comparison with results for the 77 meson production 

The comparison of the distributions of the square of the proton-proton (spp) and proton- 
meson {sp-meson) invariant masses between the pp — > pprj' and the pp — > ppr] reactions is 
presented in figure 110. 12[ For the proton-meson system the comparison was performed for 
the kinetic energy of the proton-meson system {^Sp-meson — "ip — rrimeson) and not as a 
function of Sp-meson bccausc the range of the Sp^ and Sp^' are different due to the different 
masses of the rj and rj' mesons. But the range of [y/Sp-meson — n^p — n^meson) is the same 
since the measurements for the r] and ij' production were by purpose performed at about 



Cross sections 



71 



the same excess energy. The differential distributions as a function of Spp are compared 
directly, because we compare the same system at the same energy range. 
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Figure 10.12: The comparison of the distributions of the square of the proton-proton invariant 
mass (spp) (left panel) and of proton-meson kinetic energy {^Sp-meson — rup — mmeson) (right 
panel). The distributions for the pp pprj' reaction (closed squares) were normahzed to the 
distributions for the pp — > pprj reaction (open squares) with respect to the total cross section. 



In both panels it is seen that the experimental points indicating the pp — > ppT] mea- 
surement (open squares) are in agreement with those from the pp — > pp-q' reaction (closed 
squares) within the error bars. It is unexpected that the shapes do not differ, showing the 
same enhancement at the same values of the square of the proton-proton (spp) invariant 
mass. 

The total cross sections as a function of excess energy for the pp — > ppr] and pp — > pprj' 
reactions and the comparison of the ratios of the production amplitudes for both reactions 
and the production amplitude for the pp — > ppn^ reaction showed that the interaction 
within the proton-ry should have been stronger than the proton-r/' system interaction. 
Therefore, if indeed the r^'-proton interaction is much smaller than the r/-proton as inferred 
from the excitation function, then the spectra presented here rather exclude the hypothesis 
that the enhancement (see figures and discussion in section 12.41) is due to the interaction 
of the meson and the proton. Moreover, based on those distributions it is not possible to 
disentangle between the hypothesis of the admixture of higher partial waves during the r] 
production or the energy dependence of the production amplitude. 
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10.4 Spp and Sp^/ distributions in view of theoretical predictions 

The interaction should manifest itself in the regions where particles possess small relative 
velocities as it was mentioned in previous sections. 

In figure [T0.13l the comparison of the experimentally determined cross sections as a function 
of the squared invariant mass of the proton-proton system [spp) to theoretical calculations 
is presented. One can see the strong deviation of the experimental spectrum from the 
phase space predictions (solid lines), normalized in amplitude. It can be also seen that the 
inclusion of the FSIpp [22] improves the agreement but still there is a significant discrep- 
ancy between the data and the calculationJl]. Here, it is also presented that two different 
possible normalizations of the model calculations do not agree with the shape of the exper- 
imental distribution. The left panel shows the fit to the lower values of the proton-proton 
invariant mass. The right one depicts the fitting to the high energy region of Spp. 




Spp[GeVV] Spp[GeVV] 

Figure 10.13: Distribution of the square of the proton-proton (spp) invariant mass determined 
experimentally for the pp pprj' reaction (closed squares). The integrals of the phase space 
weighted by a square of the proton-proton on-shell scattering amplitude (dotted lines)-FSIpp have 
been normalized arbitrarily to the lower values of Spp in the left panel and to the higher values 
of Spp in the right panel. The expectations under the assumption of a homogeneously populated 
phase space are shown as solid lines. 

These different normalizations are presented for the invariant mass of the proton-r]' 
subsystem in figure 110.141 



detailed description of the proton-proton FSI parameterization based on the proton-proton on-shell 
amplitude is presented in appendix [B] 
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Figure 10.14: Distribution of the square of the proton-77' (sp,,') invariant mass determined exper- 
imentally for the pp pprj' reaction (closed squares) . The integrals of the phase space weighted 
by a square of the proton-proton on-shell scattering amplitude (dotted Hnes)-FSIpp normahzed 
arbitrarily to the lower values of Spn' in the left panel and to the higher values of Spn' in the right 
panel. The expectations under the assumption of the homogeneously populated phase space are 
shown as solid lines. 



Inspecting figures 110.131 and 110.141 one can conclude that similarly to the case of the 
pp — > pprj reaction, the proton-proton on-shell interaction is not sufficient to explain the 
enhancement seen in both distributions, independent of the applied normalization. 



A better description is obtained when instead of the on-shell proton-proton amplitude 
the proton-proton FSI is parameterized by the Jost function [1,78]. 

The experimental distributions of Spp and Sp^/ are compared to the theoretical model pro- 
posed by V. Baru and collaborators [39] in figure 110.151 

The model is based on calculations in which the overall transition matrix element can 
be factorized into the primary production and the proton-proton final state interaction 
expressed as the inverse of the Jost function derived from the Bonn potential [1,39,79]. 
The production amplitude is calculated using the one boson exchange model (OBE) and 
is nearly constant in the calculations. One can see that the model normalized to the 
total cross section, hardly fit to the experimental data at the range of lower and higher 
Spp values. On the other hand, the Sp^' distribution is rather well reproduced by the model. 
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Figure 10.15: (Left) Distribution of the squared proton-proton (spp) invariant mass determined 
experimentally for the pp — > ppij' reaction at the excess energy Q = 16.4 MeV (closed squares). 
(Right) Distribution of the squared proton-77' {spr,') invariant mass determined experimentally for 
the pp ppr]' reaction at the excess energy of 16.4 MeV (closed squares). The experimental data 
are compared to the calculations under the assumption of the ^Pq —^^Sqs transition according to 
the model of V. Baru et al., as described in text and in [39,79]. 



A better agreement with the experimental data can be achieved by the implementation 
of contributions from the higher partial waves in the theoretical calculations. In figure 
110.161 the distributions of the squared masses of proton-proton (spp) and proton- r/' (spr?') 
subsystems are compared to the calculations of K. Nakayama et al. [12]. 
The calculations result from a combined analysis (based on the effective Lagrangian ap- 
proach) of the production of r] and r/' mesons in photo- and hadro-induced reactionJ§ 
[12,34,80,81]. The authors assumed in the calculations that the ^5*0 -^^ Pqs transition 
contributes beside the ^Pq -^^Sqs transition at the excess energy of 16.4 MeV. 
Interestingly, the authors claim that for the pp ppr]' reaction the dominant production 
mechanism for the NN NNif reaction is due to the excitation of the intermediate 
5ii resonance [34]. Figure [TO . 1 6 1 shows that both distributions are well reproduced by the 
theoretical calculations. 



^The general overview of the calculations proposed by K. Nakayama and group is discussed in ap- 
pendix [C] 
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Figure 10.16: (Left) Distribution of the squared proton-proton (spp) invariant mass determined 
experimentally for the pp — > ppij' reaction at the excess energy Q = 16.4 MeV (closed squares). 
(Right) Distribution of the squared proton-77' {spr,') invariant mass determined experimentally for 
the pp pprj' reaction at the excess energy of 16.4 MeV (closed squares). The experimental 
distributions are compared to the calculations based on the assumption that S and P partial waves 
can contribute as it is described by K. Nakayama et al. in [12,34,80,81]. 



On the other hand, as proposed by A. Deloff, one can explain the enhancement seen in 
the distributions by an energy dependent production amplitude, which in the calculations 
of V. Baru et al. [39,79] and of K. Nakayama et al. [12,34,80,81] was constant. 
In figure 110.171 the experimentally determined spectra of the squared proton-proton in- 
variant mass {spp) and proton- r/' (spr?') system is presented and compared to calculations 
by A. Deloff [13,82]. The calculations are based on the assumption that the production 
amplitude changes linearly with energjl. The calculations based on the standard on-shell 
approximation were modified allowing the linear energy dependence of the ^Pq — >^ Sqs 
partial wave amplitude. Other partial wave transitions were neglected in the model. 



''The phenomenological model of the differential cross section parameterization assuming the linear 
energy dependence of the production amplitude is presented in general in appendix IdI 
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Spp[GeVV] Sp^,[GeVV] 

Figure 10.17: (Left) Distribution of the squared proton-proton (spp) invariant mass determined 
experimentally for the pp — > ppij' reaction at the excess energy Q = 16.4 MeV (closed squares). 
(Right) Distribution of the squared proton-77' {spr,') invariant mass determined experimentally for 
the pp pp-q' reaction at the excess energy of 16.4 MeV (closed squares). The experimental data 
are compared to the calculations performed assuming a linear energy dependence of the production 
amplitude as proposed by A. Deloff [13,82]. 

In conclusion to this chapter one has to stress that from the above presented considera- 
tions, the determined invariant mass distributions (spp and Sp^') strongly deviate from the 
predictions based on the homogeneous population of events over the phase space. Also, 
the parameterization of the proton-proton on-shell interaction is not sufficient to explain 
the enhancement seen in the determined distributions. 

Furthermore, one can see that both calculations assuming a significant contribution of 
P-wave in the final state (K. Nakayama), or the model assuming a linear energy dependence 
of the production amplitude by A. Deloff, reproduce the data within the error bars quite 
well. Therefore, on the basis of the presented invariant mass distributions, it is impossible 
to disentangle which of the discussed models is more appropriate. 

On the other hand, it is clearly seen that the calculations assuming only '^Pq — >^ Sqs 
transition, independent of the energy, underestimate the experimental data at the higher 
values and overestimate them at the lower values of Spp as it is presented in figure 110.151 
Moreover, independently from the theoretical calculations, one can conclude from the 
comparison of Spp and y/Sp-meson presented in figure [TO. 121 that the observed enhancement 
is not caused by a proton-meson interaction, since the strength of proton-r/ and proton-r/' 
interaction is different but the enhancement in both cases is the same. 



11. Summary 



Differential distributions of the squared proton-proton (spp) and proton- r/' {spri') invariant 
masses, as well as angular distributions for the pp pprj' reaction have been evaluated at 
an excess energy of Q = 16.4 MeV. The data were measured by the COSY-11 collaboration 
at the accelerator COSY. 

The theoretical framework aiming to explain the unexpected enhancement observed in 
the invariant mass distributions determined for the pp —> ppr] reaction was presented and it 
was tried to verify one of the three challenging theories describing the observed behaviour 
by: 

• a significant role of proton-ry interaction in the final state, 

• an admixture of higher partial waves in the produced 
proton-proton-meson system, 

• an energy dependence of the production amplitude. 

On a basis of only pp ppr] data, it was impossible to judge which model could describe 
the effects. Therefore, in order to verify those hypotheses by a comparison of the invariant 
mass spectra for the production of two different mesons, the analysis of the pp — > pprj' 
reaction, presented in detail in this thesis, was performed. 

The reaction was measured using the COSY-11 detector setup designed especially for close 
to threshold measurements. The experiment was performed at a beam momentum of 
Pb = 3.260 GeV/c, which corresponds to an excess energy of Q = 16.4 MeV. 
A short description of the synchrotron COSY together with a closer look at the detectors 
of the COSY-11 apparatus used for the measurement was presented and its calibration and 
adjustment was described at length. 

In the analysed experiment two charged particles were registered. Their identification as 
protons was based on the independent evaluation of their velocities and momenta. The 
decay products of the r/' meson were not measured, therefore the number of pp — > pp-q' 
events was determined using the missing mass technique. 

For the first time ever high statistics invariant mass Spp and Sp^' spectra and angular 
distributions, determined in the close to kinematical threshold region for the pp pprj' 
reaction were presented. The distributions were corrected for the geometrical acceptance 
of the COSY-11 apparatus and were determined free from the background of the multi- 
pion production. Besides the differential distributions also the total cross section for the 
pp — > ppr]' reaction has been established, determined at an excess energy of Q = 16.4 MeV 
to be: a = (139 ±3 ±12) nb. 
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The evaluated invariant mass distributions were similar to those determined for 
the r] meson production. This result together with the shape of the excitation 
functions for the pp pprj and pp pprj' reactions, excludes the hypothesis pos- 
tulating that the observed enhancement is caused by the interaction betv^^een 
proton and meson. 

Furthermore, the determined Spp and Sp^' differential distributions do not allow to judge 
between the proposed theoretical models. Both, the hypothesis that an admixture 
of higher partial waves in the produced proton-proton-meson system and the 
proposal that an energy dependence of the production amplitude leads to the 
observed shape of the Spp and Spjj' distributions are consistent with the data at 
the achieved level of accuracy. 

The results presented in this thesis may be useful for the verification of other recently 
published theoretical models applied for the proton-proton-meson interaction, for example 
the calculations proposed by S. Ceci et al. in [83,84]. 

In order to disentangle the possible explanations for the observed differential distribu- 
tions more experimental data are required. The data base could be extended by studying 
the excess energy dependence of the shape of the Spp and Sp^' distributions, as it was done 
for the pp pprj reaction. Such studies would allow for a direct comparison of the differ- 
ential distributions behaviour for the r/' production with already determined quantities for 
the pp — > pprj reaction, in a similar manner as described in this thesis. 

A more sensitive test, whether higher partial waves are indeed important for the stud- 
ied reactions, could be gained from polarization observables. Therefore, measurement of 
analysing powers and spin correlation coefficients would be desirable for a better under- 
standing of the dynamics and interaction in the pprj and pprj' systems. 



A. Pseudoscalar mesons 



Mesons are bound states built from a quark q and an antiquark q' according to the Quantum 
Chromodynamics, abbreviated as QCD. The quarks q and q' can be of the same or a 
different flavour. Since quarks possess the spin 1/2, they can build a singlet state with the 
total spin J = or triplet states with the total spin J = 1. 

A positive parity is assigned to the quarks whereas a negative parity is assigned to the 
antiquarks. Using that convention for the bound state of qq' and denoting by L the orbital 
angular momentum of the system, the parity of the meson built out of the quark pairs is 
P= (-1)^+1. 

Prom the three lighter quarks u, d and s nine possible qq' combinations can be built. 
This allows to built a meson nonet according to the SU(3) frame including the octet and 
the singlet state: 

3®3 = 8©1. (A.l) 

The assumption that these qq' systems are the ground state combinations of the quark- 
antiquark pairs with the angular momentum L = implies, that the parity P of the con- 
structed mesons should be equal to P = — 1. With that conditions (L = and P = —1) 
mesons with the internal spin J = 1 are called vector mesons, whereas mesons with the 
internal spin J = are called pseudoscalar mesons. 

The ground states of the pseudoscalar meson nonet is presented in figure lA.ll 




L = -l I, = -1/2 I, = I, = 1/2 I, = l 



Figure A.l: The ground states of the pseudoscalar mesons (J^ = 0"). The stran ffeness S of the 
mesons is plotted versus the third component of the isospin I3. 
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The quark structure and masses of the ground states of the pseudoscalar mesons are 
listed in table lA.ll 



JT DcLlQUbL-djldji lllcoUll 


Vt^LldjiK DliU.L.lU.ic 


iVldoo iVlc V 




1 /x/2(vv — dJ) 

-L / V \ Li Li LLLL 1 


1 S4 Q8 


7r+ 


ud 


1 57 

-L <J t/ . (J J 


7r~ 


du 


139.57 


K+ 


us 


493.68 




us 


493.68 




ds 


497.61 




ds 


497.61 




A{dd + uu) + B{ss) 


547.85 




a[dd + uu) + b{ss) 


957.66 



Table A.l: List of the pseudoscalar mesons as the quark-antiquark combinations. Meson masses 
are taken from reference [5]. 



B. Parameterization of the proton-proton Final State 
Interaction 



The total transition amplitude for the pp — > ppX process can be written (in the Watson- 
Migdal approximation [85]) as [86]: 

\Mpp^ppx\^ = |Mo|2 • \Mfsi\^ ■ ISI, (B.l) 

where Mq denotes the total production amplitude, Mpsi refers to the elastic interaction 
among particles in the exit channel, and ISI stands for the reduction factor in the initial 
channel of the colliding protons. 

The element iMir^/l is predominantly due to the pp-FSI (Mpp_>pp), which can be calculated 
as [87]: 

Mpp^pp - ^— , [B.2) 

where = ^a^'^^lx is the Coulomb penetration factor [88], r/c is the relativistic Coulomb 
parameter. 

The phase shift 6pp(^So) is calculated using the Cini-Fubini-Stanghellini formula with the 
Wong-Noyes Coulomb correction [18,19,89]: 



p ctg6pp + 2p r?, h{vc) = -— + ^ - :r^^, (B.3) 

^ 1 + Qpp 



where h{r^,) = -ln{r^,) - 0.57721 + r,l Y.n=i ;i7raT [90]' 



O'pp = — 7.83 fm. 



and bpp = 2.8 fm stand for scattering length and effective range, respectively [19]. 

The parameters Ppp = 0.73 fm^, and Qpp = 3.35 fm^ are related to the shape of nuclear 

potential derived from one-pion-exchange model [19]. 
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C. Formalism of a combined analysis of photo- and 
hadro-production of the 77' meson 

Investigations of the ij and rj' meson production are performed for the following reaction 
channels: 

'J + N ^ N + M, 

IT + N ^ N + M, 

N + N^N + N + M, (C.l) 

based on a relativistic meson-exchange approach, where M = r],r]'. 
According to the model described in [91,92] the transition amplitude M [91] 
for the pp — > ppif reaction is given by: 

M = {l + TfGf)J{l + T,G,), (C.2) 

where Tjj denotes the NN T-matrix interaction in the initial (i) or final (/) state, 

and Gij is the corresponding two-nucleon propagator [91]. 

J sums all basic rj' production mechanisms depicted in figure IC.li 




P 
P 



+ 











R 


y- 

■M 



p, , p'.^y^^^^'^^Pi 

+ m:: ^ 

P2 P'2-<^ P2 

Pi p',^y''''-^—p, 

+ m;; R 

P2 P'2-<^ P2 



P r 
P'2 



ip,co 



lp,co 

m — P2 



Figure C.l: Basic production mechanisms for the pp pprj' reaction. Time proceeds from right 
to left for each individual graph. N and R denote the intermediate nucleons and resonances, 
respectively. Af stands for all possible exchanges of mesons as: tt, 77, p, u, a and ag for the nucleon 



graphs (the first two graphs of the sum in figure [C3|) and tt, p and lo for resonance graphs (the 
third and the fourth graphs of the sum in figure IC.ll ) . 
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The free parameters of the model - i.e. resonance parameters and the NNM couphng 
constant - are fixed to reproduce the available data for photo- and hadron-induced reac- 
tions listed in equation IC.ll 



D. Linear energy dependence of the production amplitude 



The following considerations and notation are adopted from reference [13]. 

The cross section for the pp ppr]' reaction can be expressed in the following manner: 

= \C{k)\^[a + b P2{p ■ q)] + Co ir])[d ReC{k) + e ImC{k)]P2{p ■ k) (D.l) 

where in the standard notation dLips stands for the invariant three body phase space ele- 
ment, p denotes the initial momentum of the proton, k stands for the relative momentum of 
the two protons in the exit channel, and q corresponds to the relative momentum between 
meson and proton-proton pair in the final state. C{k) and Co{r]) denote enhancement 
factor and Coulomb penetration factor, respectively. 

The a parameter denotes the modulus of the sole squared production amplitude which 
survives at threshold and is associated with the ^Pq -^^Sqs transition. The parameter b 
stands for the interference term between the latter amplitude and the ^Pq — ^^Sqc? ampli- 
tude. Parameters d and e correspond to the interference with a ^P2 -^^ D2S amplitude. 
The whole set of parameters depend on the final state momenta. 



Looking at the threshold behaviour, one can see that the a parameter turns to be 
constant, the b parameter is proportional to q^ and either d or e are of the order of k'^. 
In such a situation the a parameter has to be expanded to the same order by setting: 
a = oo + aiq'^ where and oi are two unknown parameters. The oq parameter can be 
absorbed in the normalization process and the ultimate expression for the cross section 
reads as: 

" oc \C{k)\''(l + -^[x + yP2{p-q)])+-^CQ{r])\zrReC{k) + z^ImC{k)]P2{p-k 



dLips V qmax ^m 

(D.2) 

where x, y, Zr and Zi are the dimensionless parameters to be determined. The x parameter 
stands for the correction of the order of q^ to the dominant transition and the y, Zr and Zi 
parameters are the measure of the D- waves admixture. 



85 



Appendix D 



Acknowledgments 

First of all, this dissertation would have not come into life without enormous help 
of Prof. Pawel Moskal. I am really honoured and have a great pleasure to work 
with Him. I admire His knowledge and I would like to express my admiration for 

His patience, lots of fruitful ideas and great concepts. 
I acknowledge Prof. Walter Oelert for being my supervisor during the stay in 

Julich. 

I am very grateful to Dr. Dieter Grzonka and Prof. Walter Oelert for all 

challenging ideas and help. 
I thank Colleagues from the COSY-11 group for nice cooperation. 
I would like to thank Dr. Rafat Czyzykiewicz for answers to all my questions. 

I really learnt a lot from Him. 
Especially, I thank Marcin Zielinski for the whole hunch of expert's suggestions. 

For careful reading and correcting this manuscript I am grateful to 
Prof. Walter Oelert, Dr. Dieter Grzonka and Prof. Pawel Moskal. 
I acknowledge a big support of my family during the time of my studies. 
Last but not least, I would like to thank my wife for being with me, 

sharing physics passion and daily life. 



References 



[1] p. Moskal, Habilitation Dissertation, 

Jagiellonian University Press (2004), e-print arXiv:hep-ph/0408162 

[2 



[3 
[4 
[5 
[6 

[7 
[8 
[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 

[17: 

[18 
[19 
[20 



S. D. Bass, Acta Phys. Slov. 56 (2005) 245. 
S. D. Bass, Acta Phys. Pol. B Proc. Supl. 2 (2009) 11. 
S. D. Bass, A. W. Thomas, Phys. Lett. B 634 (2006) 368. 
C. Amsler et al., Phys. Lett. B 667 (2008) 1. 
C. P. Jessop et al., Phys. Rev. D 58 (1998) 052002. 
G. Brandenburg et al., Phys. Rev. Lett. 75 (1995) 3804. 
K. Hagiwara et ah, Phys. Rev. D 66 (2002) 010001. 
P. Moskal et al., Phys. Rev. C 69 (2004) 025203. 

A. Fix, H. Arenhovel, Phys. Rev. C 69 (2004) 014001. 

A. Fix, H. Arenhovel, Nucl. Phys. A 697 (2002) 277. 

K. Nakayama et al., Phys. Rev. C 68 (2003) 045201. 

A. Deloff, Phys. Rev. C 69 (2004) 035206. 

S. Brauksiepe et al., Nucl. Instr. & Meth. A 376 (1996) 397. 

P. Klaja et al., AIP Conf. Proc. 796 (2005) 160. 



http:/ /ikpellOl. ikp.kfa-juelich.de/ 



R. Maier, Nucl. Instr. & Meth. A 390 (1997) 1. 
H. P. Noyes, Ann. Rev. Nucl. Si. 22 (1972) 465. 
J. P. Naisse, Nucl. Phys. A 278 (1977) 506. 
D. Sigg et al., Nucl. Phys. A 609 (1996) 269. 

89 



90 



References 



D. Sigg et al., Nucl. Phys. A 617 (1997) 526. 
P. Moskal et al., Phys. Lett. B 482 (2000) 356. 
F. Balestra et al., Phys. Lett. B 491 (2000) 29. 

R. Wurzinger et al., Phys. Lett. B 374 (1996) 283. 
P. Moskal et al., Phys. Rev. Lett. 80 (1998) 3202. 
P. Moskal et al., Phys. Lett. B 474 (2000) 416. 
A. Khoukaz at al., Eur. Phys. J. A 20 (2004) 345. 
F. Hibou et al., Phys. Lett. B 438 (1998) 41. 

E. Chiavassa et al., Phys. Lett. B 322 (1994) 270. 
H. Calen et al., Phys. Rev. Lett. 79 (1997) 2642. 

A. M. Bergdolt et al., Phys. Rev. D 48 (1993) R2969. 
J. Smyrski et al., Phys. Lett. B 474 (2000) 182. 
H. Calen et al., Phys. Lett. B 366 (1996) 39. 

K. Nakayama, Y. Oh and H. Haberzettl Acta Phys. Pol. B Proc. Supl. 2 (2009) 23. 
R. Czyzykiewicz, private communication (2008). 

P. Klaja, R. Czyzykiewicz, P. Moskal, Acta Phys. Pol. B Proc. Supl. 2 (2009) 45. 
M. Abdel-Bary et al., Eur. Phys. J. A 16 (2003) 127. 

H. Pettersson, Proceedings MENU 2007, Jiilich, Germany, 10-14 Sep 2007, 276. 
V. Baru et al., Phys. Rev. C 67 (2003) 024002. 



H.-H. Adam et al.. Proposal for WASA-at-COSY, e-print arXiv: nucl-ex/0411038 



http://wasasrv.ikp.kfa-juelich.de/WasaWiki/index.php/Main_Page 



http://www.fz-juelich.de/ikp/pisa/pisa.html 



http:/ /www. iskp.uni-bonn.de/gruppen/edda/publ/top. html 



http://users.uj.edu.pl/ ^ ufkamys/C0SY13/ 



ht t p : / / www .fz-juelich.de/ikp/ anke / de/index . sht ml 



http://www.fz-juehch.de/ikp/COSY-TOF/index_e.html 



References 



91 



[47] |http: / /www. fz-juelich.de/ikp /publications /AR2003 /en/ 



vii_collaborations.shtmn Jessica 

[48 



[49 
[50 
[51 
[52 
[53 
[54 
[55 

[56 
[5?: 

[58 
[59 
[60 
[61 
[62 
[63 
[64 
[65 
[66 
[67 

[68 
[69 



http://usparc.ihep. su/spires/find/experiments/www2?ee=C0SY-NESSI 



http:/ /www. fz-juelich.de/ikp/gem/gem. html 



http:/ /www. iskp.uni-bonn.de/gruppen/momo/momo. html 



http:/ /www. iskp.uni-bonn.de/gruppen/hires/index. html 



D. Prasuhn et al., Nucl. Instr. & Meth. A 441 (2000) 167. 
H. Stockhorst et al., IKP Annual Report (1997) 160. 

H. Stockhorst, Shriften des FZ-Jiihch, Matter and Materials 11 (2002) 176. 

R. Czyzykiewicz, Diploma Thesis, 
Jagiellonian University, Cracow (2002), IKP Jiil-4017. 

P. Moskal. Annul Report IKP-KFA (1996) 35. 

J. Przerwa, Diploma Thesis, 



Jagiellonian University, Cracow (2004), IKP Jiil-4141, e-print arXiv: hep-ex/0408016 



T. Rozek, Ph.D. Dissertation, University of Silesia (2005) IKP Jul-4184. 
R. Bilger et al., Nucl. Instr. & Meth. A 457 (2001) 64. 
P. Moskal et al., J. Phys. G 32 (2006) 629. 

P. Moskal, Shriften des FZ-Jiilich, Matter and Materials 21 (2004) 249. 
H. Dombrowski et al., Nucl. Instr. k Meth. A 386 (1997) 228. 

A. Khoukaz et al., Eur. Phys. J. D 5 (1999) 275. 

B. Gugulski et al.. Annul Report IKP-KFA(I)-1992-3. 
Bicron Corporation, Newburry, OH 44065, USA. 

G. Anton et al., Nucl. Instr. & Meth. A 310 (1991) 631. 

P. Moskal, Ph.D. Dissertation, Jagiellonian University, 
Cracow (1998) IKP, Jul-3685. 

T. Sefzick, Annual Report of IKP (2002) 14, Berichte des FZ-Jiihch, Jul-4052. 

W. R. Leo, Techniques for Nuclear and Particle Physics Experiments, 
Springer - Verlag Berlin Heidelberg, 1997. 



92 



References 



[70 
[71 

[72 
[73 

[74 
[75 
[76 
[77 
[78 

[79 
[80 
[81 
[82 
[83 
[84 
[85 
[86 
[87 
[88 
[89 
[90 
[91 
[92 



P. Moskal et al., Nucl. Instr. & Meth. A 466 (2001) 448. 

M. Zielinski, Diploma Thesis, 
Jagiellonian University, Cracow (2008), IKP Jiil-4277, e-print arXiv: hep-ex/0807.0576. 

D. Albers et al., Phys. Rev. Lett. 78 (1997) 1652. 

K. Kilian, R. Maier, S. Martin, U. Pister, (1990), 
Cooler synchrotron COSY Jiilich - User Guide. 

D. Prasuhn, private communication (2006). 

H. Stockhorst, private communication (2006). 

F. Balestra et al., Phys. Rev. Lett. 89 (2002) 092001. 

J. Ritman, Habilitation Dissertation, Giessen University (2000). 

M. L. Goldberger and K. M. Watson, "Collision Theory", John Wiley & Sons, 
New York (1964). 

V. Baru et al., Phys. Atom. Nucl. 64 (2001) 579. 

K. Nakayama, H. Haberzettl, Phys. Rev. C 69 (2004) 065212. 

K. Nakayama, Y. Oh, H. Haberzettl, e-print: arXiv: hep-ph/0803.3169. 

A. Deloff, private communication (2009). 

S. Ceci, A. Svarc, B. Zauner, Acta Phys. Pol. B Proc. Supl. 2 (2009) 157. 
S. Ceci, A. Svarc and B. Zauner, Phys. Scripta 73 (2006) 663. 
K. M. Watson, Phys. Rev. 88 (1952) 1163. 

A. Moalem et al., Nucl. Phys. A 589 (1995) 649. 

B. J. Morton et al., Phys. Rev. 169 (1968) 825. 
H. A. Bethe, Phys. Rev. 76 (1949) 38. 

H. P. Noyes, H. M. Lipinski, Phys. Rev. C 4 (1971) 995. 

J. D. Jackson, J. M. Blatt, Rev. of Mod. Phys. 22 (1950) 77. 

K. Nakayama, J. Speth and T.-S. H. Lee, Phys. Rev. C 65 (2002) 045210. 

K. Nakayama, H. F. Arellano, J. W. Durso and J. Speth, 
Phys. Rev. C 61 (1999) 024001. 



